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ABSTRACT 


The Navy is in need of an organie, inexpensive, swift method to neutralize or 
sweep waterborne mines. This thesis presents an alternative to current mine 
countermeasure technologies that fulfills this criteria—the use of the Joint Direct Attack 
Munition (JDAM) to clear a minefield. It updates the general, physics-based, six degrees 
of freedom model, STRIKE35, to predict the three-dimensional, free-fall trajectory and 
orientation of a MK-84 bomb (simulating the JDAM) through a water column. It 
accurately predicts the final detonation position relative to an underwater mine in the 
very shallow water environment. Input parameters include accurate water impact speed 
and surface impact angle of attack. Because the model results compare well with 
experimental data from the Stand-Off Assault Breaching Weapon Fuze Improvement 
(SOABWFI) Program, we analyzed the trajectory of the weapon with structural failures. 
This thesis solves for the impact speed and impact angle of attack limitations to remain 
within the Technology Transition Agreement, the detonation location for each fuze delay 
setting (to include its 20% tolerance), and the trajectory changes due to different water 
densities. This gives strike planners a tactical decision aid to clear the minefield 
accurately and efficiently with existing aircraft and weapons. 
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I. INTRODUCTION 


A, ECONOMY AND SECURITY 

The global economy in which we prosper is due solely to the free movement of 
goods and services secured by the world’s navies, most specifically, our own. The 
oceans support 90% of the world’s trade and two-thirds of its petroleum transport. In the 
joint doctrine, A Cooperative Strategy for 21st Century Seapower, the heads of the U.S. 
maritime forces accurately describe the sea-lanes and supporting shore infrastructure as 
“lifelines of the modern global economy, visible and vulnerable symbols of the modem 
distribution system that relies on free transit through increasingly urbanized littoral 
regions” (Conway et al. 2007). 

Even shutting down a single port can have devastating consequences on a nation’s 
economy. During RADM Mine’s 2005 brief to the Mine Warfare Association, he told of 
the ripple effect that closing Long Beach port caused to other sectors of the U.S. 
economy in the fall of 2002. He described the ten-day lockout as taking three and a half 
months to clear the container backlog, even after they took mitigating measures in 
anticipation of the strike. He also showed a table, which he said was “dated.” Table 1 
presents the value of cargo handled by the top ten U. S. ports, showing the enormity of 
economic damage a port closure would cause. In fact, the widely quoted estimate claims 
a shutdown of port facilities would cost the U.S. economy $1.94 billion per day (Hall 
2004). 
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Port 

$Million/day 

Los Angeles 

278.9 

Long Beach 

269.0 

New York 

221.6 

Houston 

118.9 

Seattle 

88.5 

Charleston 

86.3 

Hampton Roads 

77.8 

Oakland 

68.5 

Baltimore 

56.4 

Tacoma 

53.4 


Table 1. Top ten U. S. ports by value of cargo handled [From Hines 2005] 


While the previous examples of financial impact are specific to the U.S., they 
serve to illustrate the overall criticality of the flow of international commerce, including 
energy supplies. There is also no denying the inextricable link between a country’s 
prosperity and its security. To protect and sustain global stability and interdependent 
networks, we need to maintain control of the sea. Nevertheless, laying mines is a 
relatively simple method to interrupt or deny this control. However, ports and 
chokepoints are not the only places vulnerable to disruption. During operations other 
than peace, sea-lines of communication and beaches for amphibious landings are also 
easy targets for mining. 

B, THE THREAT 

Mines are prolific. Currently, there are over 50 countries possessing the 
capability to mine. Of these, at least 30 have demonstrated production capability and 20 
have attempted to export. Commercialization and foreign military sales programs 
currently drive the abundance of mines (Federation of American Scientists 2008). 
Because mines are a flexible and relatively simple means of denying the use of an area, 
they are the perfect “fire and forget” weapon. 

One cannot ignore a mine threat. No matter if mining is only threatened or actual 
mining takes place the same action occurs to meet the potential threat. Calculating the 
value of a minefield, one can sum up the amount of time an area was denied, the man- 
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hours invested to open the area for transit, the number of ships hit or needing repair, and 
the number of people injured or killed. As long as the effeet outweighs the cost, the 
tactic is effective. Considering the production and mine-laying cost is less than 10% that 
of removal, and it taking hundreds of times longer, mines are a very effective weapon 
(Brown 1991 and Naval Studies Board National Research Council 2001), 

When a mine explodes near a ship, the effects can be devastating and the repair 
costs exorbitant. In the past 60 years, 79% of the 19 U. S. Navy ships that have suffered 
battle damage were the result of waterborne mines (Bruhn 2006). Of these, only one was 
from a modern, more expensive mine (Naval War College 2009a). This means they are 
also attractive weapons for asymmetric warfare since they are relatively low cost - as 
little as a few hundred U.S. dollars for a simple contact mine (Zolton 1991). Figure 1 
shows the return on investment for mines in the Gulf War. Moored mines (left, grey) are 
so inexpensive that even an individual could buy them. Even the more expensive bottom 
mine (right, blue) is only a fraction of the repair price. 



Mines: Favorable Investment Return 

USS Samuel. B. Roberts (FFG-58) 


Mine Cost: $1,500 vs. Repair Cost: $96 M 


USS Princeton tCG-59t 


USS Tripoli (LPH-11 


TF*- 




$1,500 
each * 


Figure 1. Mines: favorable investment return [From Naval War College 2009a] 
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c. 


TODAY’S DEFENSE 


1, Getting on Station 


The cost of repairing ships does not solely drive this problem. The price of 
getting the correct equipment to and neutralizing the minefield can cost as much or more 
than the repair. Currently Mine Countermeasure ships (SMCM) are in San Diego, 
Sasebo, and Manama while the Mine Countermeasure aircraft (AMCM) are in Norfolk 
and Manama (Figure 2). Though these forward deployed locations are in the correct 
place strategically, limited assets obviously cannot be at every chokepoint or super port in 
the world. Since both surface and air mine countermeasures (MCMs) are slow moving 
and are near the end of their service life, they typically get transported via air or sea if an 
emergency occurs away from their current operating area. 



Sasebo 
4 MCMs 


Strait of Hormuz 


Suez Canal 


Strait of Malacca 


- I' 

San Diego k ; 
8 MCMs L-’i 


Norfolk 


HM-14 & -15 


Panama Canal 


, b Bahrain 
, ^ 2 MCMs 
I ’ 1 UM.'i/i r 


HM-14 Det 

■mi 


▲ Global Super 
^ Ports (16) 


^ World Trade 
Choke points 


Taiwan Strait 


# 


Figure 2. Chokepoints and super ports [After Naval War College 2009a] 


Military airlift is the fastest mode of transportation, but it is also the most 
expensive. The typical airframe cost per flight hour is about $5,000 for the C-17, about 
$6,000 for the C-5, and about $1,500 for the C-130 (Federation of American Scientists 
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2009). This assumes that the aireraft ean land elose enough to the area denied by mines. 
In order to aoeomplish this, we must have air superiority prior to flying an undefended 
aireraft into the area. It also means the airfield is large and strong enough for the aireraft 
to land, there is room for the offload of the MCMs, and transportation is available to get 
the SMCM to the water (Naval War College 2009e). 

The Military Sealift Command (MSC) has the quiekest eargo ships in the world, 
the Fast Sealift Ships (T-AKR). They are for rapid deployment of equipment around the 
world and are normally in a redueed operating status with the ability to get underway to 
load ports within 96 hours (Military Sealift Command 2003). Their transit speed is up to 
24 kts, whieh means it will take about 18 days to make the 8,600 nm trip from Norfolk to 
the Persian Gulf. Still, this method ean preserve serviee life and reduee transit time by up 
to 50% (USS Ardent). 

2, Once on Station 

In order for the MCMs to work an area, there must not be a threat else they need 
proteetion. They laek the ability to defend themselves, exeept for elose in armaments: 
maehine guns (A/SMCM) and grenade launeher (SMCM). If threatened by an adversary, 
MCMs are unable to eontinue their work as they are restrieted in their maneuvering. 
Though AMCMs are able to sweep for mines faster than their surfaee eounterparts, what 
they gain in speed, they give up in effeetiveness and stamina. AMCMs also require flight 
time from their base (whether a ship or an airport). Nonetheless, they are more eapable in 
the shallow water to very shallow water (VSW) environments (Figure 3) due to the 25- 
foot navigational draft restrictions of the SMCM. 
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0**p 

Watar 


Depth Regimes 


Vary-Shallow 

Watar 


Craft 

Surf Landing 

Zona Zona 



Figure 3. Undersea battlespaee and types of mines found in eaeh depth [From United 

States General Aecounting Offiee 2001] 


SMCMs sweep very effeetively in a straight line, at depths greater than 25 feet, at 
about 8-10 kn ots. The size of the area then determines how long the sweep will take. 
They are also able to hunt for mine-like-objeets, and then send an Explosive Ordinanee 
Disposal (EOD) technician or the AN/SEQ-48 Mine Neutralization Vehicle (MNV) to 
inspect or destroy the object (Gasperini 2003). This shows that time is the biggest factor 
in how well an area can be cleared. The measure of success is then a probability, e.g. a 
95% confidence that the area has a 67% clearance rating. There is always a balance 
between time and risk, so a constrained timeline causes the acceptance of more risk. 

Many options exist to neutralize mines, but all options have advantages as well as 
disadvantages. Eor example, SMCMs are extremely effective but are slow and have draft 
limitations that do not allow them to operate in shallow water. SMCMs have the 
AN/SLQ-48 Mine Neutralization System (an autonomous submersible for moored and 
bottom mines), but these are also slow. AMCMs can tow the very capable sled (with the 
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MK-103 through MK-106 installed) into shallow water, but are unable to work in low 
visibility or at night. EOD teehnieians are exeellent in the VSW, but their limitations 
eome from fatigue, water temperatures, and water depth. Beyond the risk they are taking 
being in the water, if the water is murky enough to restrict vision, their risk increases. 
Marine mammals, though excellent at hunting mines and currently our only asset for 
detection of buried mines, do not neutralize mines due to the risk involved with handling 
explosives (Bruhn 2007 and Naval Studies Board National Research Council 2001). In 
order to reduce the risk to personnel and to reduce the sweep timeline without sacrificing 
effectiveness, a balance of assets is currently the best available option (Naval War 
College 2009a). 

D, TOMORROW’S DEFENSE 

In order to further decrease MCM response time and increase MCM capabilities, 
the Navy is currently integrating MCM forces, both physically and doctrinally, into all 
Navy Joint Task Forces. This reduces the reliance on dedicated MCM forces in the early 
stages of a conflict (Mullen 2000). Future dedicated capabilities include the Fittoral 
Combat Ship (FCS) to replace the Avenger class (SMCM) and the MH-60S Knighthawk 
to replace the MH-53F Sea Dragon (AMCM), just to name a few (Hines 2005). These 
systems, among others, had an initial forecast completion of FYIO, but according to 
RADM Hines (2005), the “initial program introductions continue to slip to the right.” 

The future of mine warfare (MIW) is at the mercy of technology. What happens 
to our fleet of SMCM and AMCM if an “unsweepable” mine is developed? Currently, 
modem mines have delayed arming devices, ship counters, clock-controlled live periods, 
and the ability to classify mine-hunting ships (Hartmann and Tmver 1991). This paper 
proposes adding an off-the-shelf asset to the mix of current mine neutralizers that is not 
effected by the limitations of the “unsweepable” mine. The capability of the Joint Direct 
Attack Munition (JDAM) Assault Breaching System (JABS), readily available in any 
Strike Group, is an essential addition to the MCM solution. 
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II. RESEARCH OBJECTIVE 


A, PREVIOUS 

Because of a lack of immediately available mine counter measures in any 
large quantities, the Navy is evaluating the JDAM Assault Breaching 
System (JABS) for its capability in countermine and assault breaching. 

The JDAM accuracy, repeatability (guaranteed delivery patterns for 
multiple deliveries) and fuzing options make the JABS a prime contender 
for an interim capability. Combined with bomber range and payload 
capability, this weapon system vastly improves ioint operations. 
(Almquist 2008) 

To bridge the gap between developing and aging systems, the Office of Naval 
Research (ONR) tested the effectiveness of JABS as a mine neutralizer in the surf and 
beach zone. JABS transitioned to PMS-495 and the Navy fielded it as a non-material 
solution to an existing naval challenge. This method seems to be the ultimate in naval 
organic weaponry because every carrier and amphibious assault ship’s strike aircraft can 
carry this weapon. The goal of current research at ONR is to extend the capability 
demonstrated in the surf and beach zones to the very shallow water (VSW). The new 
ONR program, called the Stand-Off Assault Breaching Weapon Fuze Improvement 
(SOABWFI) program, is testing the viability of JABS in water depths from 3.05 meters 
to 12.19 meters. 

1. JABS 

During the initial program, ONR concluded that JABS was competent in mine 
and obstacle clearance on the beach and in the surf zone (0 to 3.05 m). With accurate 
targeting information, it is effective, day or night, against unburied mines and obstacles in 
the surf and beach zones. Ultimately, it takes the man out of the minefield during 
neutralization (Almquist and McLaughlin 2008). Weaknesses include limited lethality 
against buried mines in the surf and beach zones, and a tail section that is more weakly 
constructed than the warhead resulting in tail fins or the whole tail section breaking off 
during initial water entry or during tail slap with the water cavity during its descent 
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through the water eolumn (Almquist and MeLaughlin 2008; Gefken 2006). JABS 
trajectory analysis showed that tail sensitivity is dependent on the impact angle of attack 
(AOA), defined here as the angle between the body axis and the velocity vector. If 
impact AOA is less than two degrees, the belief is that the tail will maintain structurally 
integrity beyond 28 msec. Above this limit, structural failure is likely (Watson et al. 
2009). 


2. SOABWFI 

The SOABWFI program assesses the trajectory and the lethality of JABS in the 
VSW. Again, JABS consists of the 2,000-pound MK-84 general munition with the 
JDAM kit. The GBU-31 is the weapon’s designation (Figure 4). It also gathers data on 
fuzing options and detonation times to analyze the weapon’s lethality against various 
mines. If ONR is able to prove this concept, the Navy will have a truly organic capability 
to clear VSW mines with the standoff of up to about 15 nautical miles (Naval Air 
Weapons Center, Weapons Division 2002). 



Figure 4. GBU-31(V)2 with MXU-735 nose cone [From NAVAIR 2009] 

a. Evaluation of Precision-guided Bomb Trajectory through Water 
using Scale-Model Experiments 

ONR contracted SRI International of Menlo Park, California to perform 
19 scale model experiments to determine the trajectory of the MK-84 through the VSW. 
They used this data to predict maximum stable penetration depth, to gain insight on 
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required fuze delay time for optimum target lethality, and to estimate the final detonation 
position relative to the target position. The SRI report states that past researeh shows the 
water trajeetory varies based on “initial velocity . . . and entry angle,” so they looked 
closely at these specific parameters (Gefken 2006). 

In order to examine correctly the MK-84 trajectory, they created four 
different configurations of the scale model to represent possible damage levels that may 
occur during water entry or because of tail slap within the cavitated region. The different 
configurations included a complete warhead section with (A) a tail section and four fins, 
(B) a tail section and two perpendicular fins, (C) a tail section and no fins, and (D) no tail 
section (Figure 5) (Gefken 2006). 


Figure 5. 



AJ 7t)7$>T6 BCidy and Tai 0n>CopperPiu 0 0.034i-bffi’EpoxyPiiler pmoia. p 2 .ii In. 
(A) UKS4 WARH&kO WTH TAL SECTION AND FOUR RNS 



024S4> Capper Plug 

(B) UW WA»€AO W!H TAL SECnCN AND TWO F«S 




0.2e0« Copper nug 

KO HCW Wi^HEAO WTH TAL SeOTION AND NO FNS 




{Q U<»t VW&CAO wmi NO TAL S80TKN 


l/12th scale MK-84 model design details [From Gefken 2006] 


What SRI found during its experiment was that the tail section and fins 
play an important role in stabilizing the underwater trajectory (Figure 6). The more tail 
fins the bomb lost, the greater the increase in horizontal nose position, i.e., horizontal 
displacement from the point of water entry. When the weapon lost its tail altogether, it 
tumbled soon after water entry. They also noticed the weapons with a tail section had 
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similar overall trajectory shapes (single direction curve). Their results also showed that 
the horizontal displacement increased with increased impact velocity (Gefken 2006). 



Figure 6. Trajectories for MK-84 warhead with different tail configurations [From 

Gefken 2006] 

SRI confirmed that the bomb generated and traveled in a column of air 
(called cavitation) due to the high water-entry velocity during a majority of its 
underwater travel. In all 19 of their experiments, the cavitation generated in the water 
column consisted of: 

... a tapered cone that has a maximum full-scale diameter at the end of 
the tail section of about 3.0 ft. This is about a factor of two larger than the 
maximum bomb diameter of 1.5 ft. This cavity shape was about the same 
for all of the initial water-entry velocities [and bomb configurations] 
between 392 ft/s and 991 ft/s (Gefken 2006). 

In this cavitated column, the pressure drag (produced by turbulent low-drag conditions 
down to a velocity of about 7.8 m/s for the scale model) dominates the dynamics of the 
bomb. Below this velocity, laminar flow dominates causing viscous forces on the body 
(Gefken 2006). 

Due to the terminal velocity prior to detonation being greater than 7.8 m/s, 
the SRI report looked specifically at pressure drag. They found that the drag coefficient 
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(Cd) will vary depending on whether the trajectory is linear or curved. In the linear case, 
is typically less than 0.1 for all tail configurations. For the curved trajectory, it is 
mostly between 0.2 to 0.3 with the bomb intact, but increases to 0.6 to 1.0 with two fins 
or tail without fins. Their data lastly indicates no correlation between a deviation of 
several degrees from a 90° water-entry angle and its ultimate trajectory (Gefken 2006). 

The Reynolds number (Re) for the 1/12*’^ scale must equal the MK-84 Re 
in order for the scale model to represent an accurate depiction of the trajectory of the 
bomb. In this case, the Re’s are not equal because of the inability to reduce the medium’s 
kinematic viscosity by the scale factor. Instead, because the Cd tended to be very small 
in the region of turbulent flow, the assumption made was that the model provided a good 
simulation of the MK-84 for velocities above 7.8 m/s (Gefken 2006). 

b. Underwater Bomb Trajectory Prediction for Stand-off Assault 
Breaching Weapon Fuze Improvement (SOABWFI) 

The Naval Postgraduate School used the SRI data to create a 6 degree of 
freedom (6-DOF) computational bomb maneuvering model called STRIKE35. It predicts 
the trajectory, velocity, and orientation of a MK-84 inside the water column. In order to 
accomplish this, they found the diagnostic relationships using three momentum equations 
and three moment of momentum equations, and the relationship between the drag (C^), 
lift (C/), and torque (Cm) coefficients for the rigid body position and orientation. The 
basic dynamics of the model follows. 

(1) Position and Orientation. The coordinate system, with the 
unit vectors (i, j) in the horizontal plane and (k) in the vertical direction, is earth-fixed. 
The bomb falling through the water column is assumed to be an axially symmetric rigid 
body where the head and tail coordinates are represented by r/,(t) and rft), respectively. 
The unit vector representing the body’s main axis direction (Figure 7) is e, where 
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y 

X 

Figure 7. Position vectors (r/j, r,) and the body axis unit vector (e) [From Chu et al. 

2010 ] 

In Figure 8, the center of mass (om) and center of volume (ov) are 
located on the main axis, and the distance between them (cr) is positive (negative) when 
the direction from Ov to Om is the same (opposite) as the direction of e. The lifting force 
and drag force (discussed in the next section) are F/ and F^, respectively. The angle of 
attack (a) is the angle between e and the unit vector of the rigid-body velocity (e^) 

a = cos ' (e„ •e) . (1) 

The angle between e and the horizontal plane is fi 

P = sin'(e»k), (2) 

and Y is the angle between e« and the horizontal plane 

y = sin^'(e„*k). (3) 

In STRIKE35, the location of the body is represented by the position of Om, 

r(t) = xi +yj + zk. (4) 

The body velocity (also called translation) is then 



14 



dr 

—^ = u, u = i7e 

dt 


( 5 ) 


where {U, ) are the speed and unit vector of the rigid-body velocity. Figure 9 shows the 

unit vector e„ broken into its components, e„, e^, e,7 where is the projection of 

e„ on the xy-plane 

e„ = cosy cos +cosy sin ^j +sin yk (6) 

where y/ is the azimuth and y is the elevation angle. The unit vectors are also 
perpendicular to one another, T -L e^, T , 





Figure 8. Angles, centers and forces of the bomb [Chu et al. 2010] 
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Figure 9. 



(2) Lift and Drag Coefficients. In order to define the lifting 
force, F/, they identify the unit vector () as being perpendicular to both e and e„, 


F/ is in the plane constructed by the vectors e and e„, and is in the same direction 
ase|^ X e„ (i.e. perpendicular to both unit vectors). Fj is in the opposite direction of u. 

Both forces are exerted at Oy, and are described by 

K = -fd^u > F/ = //C,, e, = el) X e„, (8) 

where fd and// are the magnitudes of the forces. The lifting force is perpendicular to the 
drag force, and their magnitudes are represented by the drag law, 

L=\c,PKU\ f,=\c,pA,p\ (9) 

Here, p is the density of water and Ayy is the underwater projection area. Cd and Ci are 
the drag and lift coefficients that are determined by substituting the SRI position and 
orientation data into these theoretical diagnostic relationships and working backwards to 
statistically solve the semi-empirical formulas of the coefficients. Lastly, U is the water- 
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to-body relative speed. Because the water velocity (Vw) is two orders of magnitude 
smaller than the bomb’s velocity (u), 

V = V„-u«-u = -t/e„. (10) 


The bomb is a fast moving, rigid body with four fins that also have 
lift and drag components. In order to solve for these values simplistically, NFS assumes 
the fins are NACA0015 airfoils, which have known drag and lift coefficients. The 
additive effect of the combined drag and lift forces on a pair of fins (two fins in the same 
plane) is F/ , so the total forces of the bomb’s fins is uF/ , where n is the number of fin 
pairs. 

(3) Momentum Equation. Differentiating equation (5) with 
respect to time gives the acceleration of the rigid body. 


du dU rrde„ 

- = -- ■ 

dt dt dt 

This makes the momentum equation in the earth-fixed coordinate system 


m 


dU 


de,. 


^u+U- 

V dt dt 


where F is the resultant force consisting of 

F = F + F, + F. + F + uF^ . 

g b a I c 

The gravity and buoyancy forces are 


( 11 ) 


( 12 ) 


(14) 


F^ = -mgk, F^ = pDgk, (15) 

where m is the bomb mass, g is gravity, and fl is the volume of the rigid body. 

Computing the inner products between equation (12) and the unit 
vectors () for n = 0 leads to 


= m)gk.e„ - , 

dt 


(16) 
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( 17 ) 


mU^ = [{pU-m)g\^ + 


and 


mU cosf^^ = [(pn-m)gk + /,e,]*e^. 
dt 


(18) 


The vector [(plT - m)gk + f^Q^ ] is 


[(pn-m)gk + /,e,] = (m^ + /Je„ +ml7^e^+mt7 cosy^el- (19) 

dt dt dt 


dU 


dy • 


Because m-h/^ is in the e«-direction, mU — is in the -direction, and 


dt 


dt 


mU cos/^^ is in the -direction, they used the perpendicularity condition and the 
dt 


inner product of (19) by the unit vector C/ to calculate 


[(-m + pn)^k + /,e,]*e, = 


mil— el +mU cosy^^e^ 
dt “ dt " 


•e,- 


( 20 ) 


(4) Moment of Momentum Equation. The moment of 
momentum equation, relative to Om, is 


^ dO. 

J*—= M,+M,. 
dt 


( 21 ) 


Here J is the gyration tensor, which is a diagonal matrix in the body-fixed coordinate 
system. 


J = 


J, 0 0 

0 0 

0 0 /, 


( 22 ) 


J\ , J 2 , and /s are the moments of inertial, and for an axially symmetric, rigid body J 2 
and J 3 are equal. The rotation velocity of the rigid body is 


^=f!5^e+n,(nx.)+^, 

dt dt dt 


(23) 


where n*is the rigid body’s angular velocity vector, and its components are the bank 
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angle along the unit veetor e, and the other azimuthal and elevation angles 

perpendicular to e ( H ). The torque by the buoyancy force (or buoyancy torque) is 

= -ere X (yorigk) . (24) 

Lastly, the hydrodynamic torque (M/,), relative to Om, is computed 

M, = -ae X (F^ + F,) - X F/ + . (25) 

Here is the anti-translation torque by the moment of drag and lift forces and is 
the antirotation torque, which is dependent on the Re. With a little manipulation and the 
assumption that this rigid-body has very slow or no self spinning (Q^^^O), the 
component of the moment of momentum equation that is perpendicular to e is 

Jj—— = Ms-cr/?ngexk-o-(/^exe„+/,exe,)-no-^.//exef +M,^+M^. (26) 

at 

(5) Hydrodynamic Coefficients. Ultimately, they rewrote the 

coefficients as 


2[(pn-m)gk*e„ 


■m- 


C,= 


dU. 
dt ' 




(27) 


dy 


dy/ 


2 [mU(el — + el cos / —) • e, - (pU - m)gk • e, ] 


Q = 


dt 


dt 


pAP^ 


(28) 


and 


= 


^2 ^ • el + erpUg (e X k) • e' 
dt _ 


+ Cf(nUt^cos;'-^^)e,.e‘+-^[C^(exe.)»e*+C,(exe,).e‘] 
K U dt L 


(29) (Chuet al. 2010). 
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(6) STRIKE35 Verification with 1/12* Scale MK-84 Models. 
STRIKE35 uses the Boeing Corporation’s specifications data of the full-scale MK-84 
(length, mass, mass inertia properties) for all variants (four, two, and no fin, as well as no 
tail). It solves the momentum equations and the moment of momentum equation using 
the numerical method and the derived hydrodynamic coefficients. Successful verification 
occurred by inputting the SRI configuration, water impact speed, and impact angle data 
to compare the STRIKE35 predicted trajectory and orientation to the observed SRI runs 
(Figures 10 through 13) (Chu et al. 2010). 


(a) (b) (c) (d) (e) (f) (g) 



X(m) X(m) X(m) X(m) X(m) X(m) X(m) 


Figure 10. Trajectory comparison of predicted and observed MK-84 with tail and four 
fins with initial water entry conditions: (a) 132 m/s and 0.8° AO A, (b) 297 m/s 
and 0° AO A, (c) 295 m/s and 1.2° AO A, (d) 302 m/s and 1.5° AO A, (e) 227 m/s 
and 2° AOA, (f) 219 m/s and 1° AOA, and (g) 119 m/s and 1.8° AOA [From Chu 

et al. 2010] 
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(a) Initial U = 295m/s 


(b) Initial U = 290m/s 


(c) Initial U = 297m/s 



Figure 11. Trajectory comparison of predicted and observed MK-84 warhead with two 
fins and initial water entry conditions; (a) 295 m/s and 0° AO A, (b) 290 m/s and 
0° AO A, and (c) 297 m/s and 0° AOA [From Chu et al. 2010] 


(a) Initial U = SCWm/s (b) Initial U = 299nrVs (c) Initial U = 291 m/s 



Figure 12. Trajectory comparison of predicted and observed MK-84 warhead with tail 
section and no fin with initial water entry conditions; (a) 304 m/s and 0° AOA, 
(b) 299 m/s and 3° AOA, and (c) 291 m/s and 1.9° AOA [From Chu et al. 2010] 
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(a) Initial U = 296m/s 


Expenment 3D Model 



Figure 13. Trajectory comparison of predicted and observed MK-84 warhead with no tail 
section and initial water entry conditions: (a) 296 m/s and 4.3° AOA, (b) 301 m/s 
and 0° AOA, and (c) 300 m/s and 1.3° AOA [From Chu et al. 2010] 

The paper on the creation of STRIKE35 concludes that Cd and C/ depend 
less on fl and more on Re and AOA. They found that Cd increases with Re, for constant 
AOA. Holding Re constant, Cd monotonically increases with AOA from vertical to a 
horizontal drop (0 to 90° AOA). The depends only on Re and AOA. C/ also depends 
on Re and AOA, but the reliance is more complicated (Chu et al. 2010). 
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B, CURRENT ANALYSIS 

The creation and verification of STRIKE35 used the same data. As seen in 
Figures 10-13, the trajectory of the scaled down MK-84 is represented well using the 
computer model. This thesis amends STRIKE35 in order to validate the program’s 
physics for live and inert GBU-31s dropped into a 7.62 m and 12.2 m test ponds. With 
the newly verified model, we solve for the limiting impact AOA and impact speed to 
remain within a Technology Transfer Agreement (TTA) requirement. Last, this thesis 
finds the limiting impact AOA and impact speed for all the delay settings, to include the 
20% tolerance in the weapon specification document, in the FMU-139 fuze in an ocean 
environment (Boeing 2009). 

By knowing the underwater trajectory of a GBU-31, the Navy can usefully 
employ this truly organic, proven warhead efficiently and effectively to quickly clear a 
minefield without harming people or trained marine mammals. The JDAM is a joint 
asset that can neutralize a mine any time of day and in any weather condition. This thesis 
provides an accurate, cost effective tool for determining the final detonation position 
relative to the water impact location. 
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III. SOABWFI FLIGHT TESTING 


A. TEST PREPARATION AND PROCEDURES 

1. Test Ponds and Targets 

Naval Air Warfare Center, Weapons Division (NAWC/WD) ereated two eireular 
ponds in the middle of Indian Wells Valley, California (Figure 14). They are built the 
same, both having a eireular bottom with a diameter of approximately 30.5 m, exeept for 
the depths. The smaller pond is about 7.6 m deep and the larger is about 12.2 m deep. 
Sloping sides (2:1) ereate a surface diameter of roughly 61 m for the smaller pond and 79 
m for the larger diameter. Both ponds have a ramp built into the side for vehicle access. 
A plastic liner covers the dirt to contain the brackish water (Figure 15) that is supplied by 
a 206 m deep, on-site well that fdled both ponds at about 800 gallons per minute (Figure 
16) (Watson et al. 2009). Placed inside the water are fully operational, moored, foreign 
mines filled with simulant, instead of TNT. Their nomenclature is “JV” as to avoid 
classification issues. 



Figure 14. Artificial ponds used for the flight tests at NAWC/WD [From Watson et al. 

2009] 
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Figure 15. (a) Man-made lake (approximately 7.6 m deep) eovered in plastie with the 

side road (b) and four JV mines [From Thomas Kelsh 2009] 



Figure 16. Ground water from a PVC pipe fills the test pond 

2. Instrumentation 

These experiments use two types of high-speed digital eameras, light sensors, 
pressure sensors, and a global positioning system (GPS) to colleet the data. The range 
eameras eapture 60 frames per seeond and the two Phantom eameras eapture 1,000 
frames per seeond. From the shore of the pond, they eolleet data that assists in solving 
for the loeation, speed, and orientation of the weapon at the time of water impaet (Figure 
17). Using orthogonal images from the Phantom eameras, they eould also solve for the 
water impaet AOA. The light and pressure sensors provide the time and depth of 
detonation for the inert weapons equipped with a fuze and booster. The booster fires at 
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the same time as the fuze, sending out a pressure pulse and light flash that is pieked up by 
the sensors. This, of eourse, assumes that the speed of light is infinite and pressure 
translates short distanees almost instantaneously in ineompressible fluid (Boeing 2009a). 
The horizontal displaeement of the weapon in the pond is determined by eomparing the 
distanee between the water impaet and pond bottom impaet. The eamera images help 
solve for the water loeation and the Trimble 5800 GPS system loeates the pond bottom 
impaet loeation surveying the holes (Figure 18) (Watson et al. 2009). 



Figure 17. High speed digital eamera 



Though the image quality of the Phantom eameras is superior to the range 


eameras, their orientation aeouraeies are unknown due to their non-surveyed plaeement. 
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Preliminary findings from Underwater Trajectory Test Three (UTT-3) state that the 
vertical axis alignment is good to only about 2°, making it unlikely to find absolute 
measurements of body axis angle or velocity vector direction at impact. Absolute errors 
in the camera measurements cancel out when looking at the differences between the two, 
though. Therefore, the AOA (defined here as the difference between the body axis angle 
and the velocity vector) is the accurate comparison between the individual analyzed runs 
(Boeing 2009c). 

3. Aircraft and Weapons 

An F/A-18F Super Hornet, proceeding at 0.80 Mach, dropped live and inert GBU- 
31s from 35,000 feet mean sea level (MSL). Release occurred approximately 5-7 miles 
from the pond in order to give the glide weapon enough kinematic energy to orient itself 
vertically above the designated point of impact (DPI). The desire is to have the velocity 
vector aligned with the munitions axis (zero AOA), and both vectors perpendicular to the 
fiat, water surface. All of the GBU-31s penetrated the water within the prescribed 
delivery error of less than 2 m Circular Error Probable (CEP) at velocities between 382.5 
and 394.9 m/s. 

Every JABS in the experiments had the MXU-735 nose cone and the tail 
telemetry (TM) kit installed. The nose cone is blunter than the Navy’s Ogive nose used 
in previous experiments (Eigure 19). The flatness forces a larger cavitation tunnel for the 
weapon to proceed through. The TM provides data, via line of sight transmission, on 
various flight parameters such a velocity, heading, altitude, and angle of attack. Since 
there is not a line of sight from the pond to range control, the TM’s lowest data 
transmission was about 32.9 m above the pond. 


28 



MXU-735 


Na\’y Ogive 

Figure 19. Differing nose cones between JABS and SOABWFI [From Boeing 2009c] 

The weapons that had fuzes were equipped with an FMU-139 B/B with available 
delay settings of 0, 10, 25, and 60 msec. Selection of the delay depends on which types 
of targets the weapon is to attack. The explosive in the live weapon is PBXN-109, 
whereas the inert weapons have filling to maintain appropriate weight and balance. 

B. UNDERWATER TRAJECTORY TEST ONE (UTT-1) 

The following description includes all that is not standard about the first 
Underwater Trajectory Test of the SOABWFI program. The primary purpose of this 
experiment was to collect data to measure the underwater trajectory deviation of the 
JABS traveling through 7.62 m of water. A secondary objective was to obtain 
information to help characterize the error budget in fuze delay, which is nominally ±20% 
according to the weapon specification document (Boeing 2009b). 

On August 7, 2008, the Super Hornet dropped four inert GBU-31s into the 
shallow pond. The delivery error, defined as the (x, y) distance between the DPI and the 
actual weapon impact point, averaged 1.3 m in this test. Three of the four weapons were 
equipped with a live FMU-139 fuze (and booster) set for a 10 msec delay, though the 
sensors internal to the pond failed to detect the flash of light and pressure pulse caused by 
the fuze booster (Boeing 2009b). 






A repeater was not available for retransmission of the TM signal to range control, 
so the last receipt of the line of sight signal was at best 32.9 m above the pond. TM data 
was able to confirm that the mission planning for the target elevation was incorrectly set 
for 7.62 m below the bottom of the pond instead of the desired 7.62 m above (at the water 
surface). This human error likely introduced some uncorrected weapon AO A at water 
impact; this is because the weapon sets the programmed AOA directly before impact, but 
the weapon also thought it had 15.24 m of descent left to set the AOA (Boeing 2009b). 

Table 2 shows that only two of the bombs have complete water impact 
information. For the first drop, the last known velocity was at about 2.0 km above 
ground level. The other three drops, in order, had TM data for 32.9, 55.5, and 78.0 m 
above the pond (Boeing 2009a). The fourth drop had one of the Phantom cameras fail. 
Typically solving for the AOA takes two orthogonal vectors. We know the total AOA 
must be greater than the largest component and smaller than the sum of the components. 
From this, we know that having only one vector we can only solve for the lower bound 
(Boeing 2009b). Although they made no effort to find and catalog debris from the 
weapons, there was some evidence that at least two fins came off (Boeing 2009b). 


Run 

Impact 

Velocity 

(m/sec) 

Angle of Attack 
(degrees from horizontal, 
±0.5 degrees) 

Total Nose 
Displacement 
(m from surface impact) 

Final Nose 
Depth 
(m) 

1 

unk 

1.66 

0.11 

25 

2 

389.4 

1.44 

0.66 

25 

3 

388.2 

1.59 

0.56 

25 

4 

388.8 

>1.91 

0.57 

25 


Table 2. UTT-1 known characteristics 


C. LIVE FIRE FLIGHT DEMONSTRATION ONE (LFFD-1) 

On September 9, 2008, the first live fire demonstration of the JABS against mine¬ 


like (JV) targets occurred at the shallow pond. Because of an incorrect fuze delay setting 
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(human error; zero was set), the bomb detonated prematurely. The test results were 
invalid because the weapon detonated very shortly after impact. A duplicate assessment 
(LFFD-1 Repeat) then occurred on November 19, 2008 with satisfactory results; the 
description below is all that is non-standard about the repeat test (Boeing 2009a). 

The primary purpose of the experiment was to demonstrate the lethality of one 
GBU-31 against several mine threats normally found at 7.62 m water depth, which are 
generally different from those previously tested during JABS. A secondary objective was 
to obtain pressure measurements in order to measure the explosive output and correlated 
damage to the targets. Since the pressures and correlating damage are incidental to this 
thesis, as well as classified, the forthcoming analysis will not discuss them (Boeing 
2009a). 

The circle in Figure 20 represents the flat bottom of the small pond with the DPI 
at the center (0, 0). Plotted on this figure are the locations of the LFFD-1 test mines with 
East in the positive x direction and North in the positive y direction (m). Two mines, 
located at (0, ±7.3), had a maximum predicted neutralization distance of 15.2 to 18.3 m 
by the JABS. Two mines, located at (±9.8, 0), had a maximum predicted neutralization 
distance of over 30.5 m by the JABS. The remaining three mines, at standoff distances of 
3.0, 4.6, and 6.1 m from the center (2.1, -2.3), (3.5, 3.1), and (-5.6, 2.2), are typically 
found at 4.6 m of water and have a maximum neutralization distance of 6.1 to 9.1 m by 
the JABS at this lesser depth (Boeing 2009a). 
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Figure 20. Florizontal positions of JV mines in the pond [From Boeing 2009a] 

Analyzing the NAWC/WD eamera data, the weapon hit the pond within a few 
meters (1.3, -0.6) of its DPI (0.0, 0.0). Unlike the first LFFD-1, this weapon detonated 
approximately 12 msee after impaet (Figure 21). Post explosion, the lake had lost 
approximately 0.9 m of water, and a grey film sat on the perimeter (Figure 22). Large 
pieees of debris were outside of the pond, some as far as 30.5 m away from the edge 
(Figure 23). Visual inspeetion of the JV mines showed damage, and draining the pond 
(Figure 24) revealed a large, bottom erater, presumably eaused by the effeet of the bubble 
expansion seen in Figure 21, frames three through six (Boeing 2009a). 
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Figure 22. Post detonation test pond 
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Figure 23. Test pond debris 



Figure 24. Drained test pond 
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Table 3 shows the repeat test’s pertinent water impaet eharaeteristics. Using a 
software program ealled LOCATE, the Naval Surfaee Warfare Center, Carderoek 
(NSWC/CD) solved numerieal equations for the weapon position and time at detonation. 
An assumption made to ealeulate this position was that the eenter of the eharge was 
approximately at the eenter of the 2.4 m long explosive seetion, which is about 1.2 m 
from the tip of nose. Results for the position at detonation are (1.5, -0.6, -3.1) m. The 
detonation x- and y-coordinates are very close to the impact position (1.3, -0.6, 0.0) m. 
Based on the digital imagery, the actual time of detonation was 12 msec (Boeing 2009a). 
Using the Pythagorean Theorem to calculate the actual travel distance of 3.6 m, we find 
the detonation time and distance correlates well; (382.5 m/sec)(0.012 sec) = 4.59 m 
(nose) - 1.2 m (explosive) = 3.39 m. This does not account for the deceleration of the 
weapon penetrating the water, but allows the reader to see that the assumptions made still 
compute to a reasonable error of 0.21 m. 


Run 

Impact 

Velocity 

(m/sec) 

Angle of Attack 
(degrees from 
horizontal, ± 0.5 
degrees) 

Total Explosive 
Displacement 
(m from surface 
impact) 

Final 
Explosive 
Depth (m) 

Time to 
Detonation 
(sec) 

1 

382.5 

2.04 

0.18 

3.1 

0.012 


Table 3. LFFD-1 known characteristics 


D. UNDERWATER TRAJECTORY TEST TWO (UTT-2) 

The subsequent explanation is the non-standard occurrences from the second 
Underwater Trajectory Test of the SOABWFI program, which took place in the deep 
pond. The primary purpose of this experiment was to collect data to measure the 
underwater trajectory deviation of the JABS traveling through 12.2 m of water. A 
secondary objective was to evaluate the fuze function with a 25 msec delay, and collect 
pressure data during the underwater trajectory (Boeing 2009c). 

On February 28, 2009, the Super Hornet dropped two inert GBU-31s into the 
deep pond. The original plan was to drop four, but the third weapon hung on the bomb 
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rack, requiring the pilot to exeeute hung ordnanee proeedure. This ultimately terminated 
the flight without dropping the last two bombs. A strong wind event had oeeurred the 
previous day, so the depth of the pond was about 11.9 m. During the two drops, the wind 
averaged 1.5 m/s, whieh was too low to have had a major impaet on the weapons’ AOA. 
The average delivery error of the two drops was 1.5 m. As well, the seeond weapon had 
a live FMU-139 fuze (and booster) set for a 25 msee delay (Boeing 2009b). 

This experiment ineluded four towers eaeh with four underwater video eameras, 
whieh reeorded (at 60 frames per seeond) the first drop aimed at the West aim point (AP) 
(Figure 25). These video images elearly show this weapon did not eontinue straight 
beyond about 6 m depth. In faet, they show evidenee that the tail of the weapon hit the 
eavity wall when the nose was about 5.5 m depth and broke off at least a pieee of it 
(Figure 26). It is elear that the interior of the tail seetion was fully exposed to the water, 
but unelear from the photos if the tail separated fully from the body before bottom impaet 
(Boeing 2009b). 



-50 40 -30 -20 -10 0 10 20 30 40 50 


Figure 25. Aim points (AP) and eamera towers [From Boeing 2009b] 
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Figure 26. Tower A, camera 2 images showing tail parts leaving the cavity (sequential 

photos with 66 msec between) [From Boeing 2009b] 

Table 4 shows the second underwater trajectory test’s pertinent water impact 
characteristics. Both weapons had such an angle when hitting the pond’s bottom that it 
was difficult to locate where the nose first penetrated the pond bottom. To reconstruct 
where the hole was located in the plastic liner, to within 0.3 m accuracy, they had the 
liner extracted and stretched out. The second weapon, aimed at the North AP, hit the 
bottom with shallow enough angle to burrow under the liner, climb the North pond face 
for a bit before re-entering the water going upwards, and finally getting airborne (Figure 
27). 



Figure 27. UTT-2(2) get airborne after penetrating the pond bottom [From Boeing 

2009b] 
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Run 

Impact 

Velocity 

(m/sec) 

Angle of Attack 
(degrees from 
horizontal, ±0.5 
degrees) 

Total Nose 
Displacement 
(m from surface impact, 

± 0.3 m) 

Final Nose 
Depth 
(m) 

1 

392.2 

1.27 

0.72 

39.1 

2 

394.9 

2.27 

0.43 

39.1 


Table 4. UTT-2 known characteristics 
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IV. MODEL VERIFICATION USING JDAM AND APPLICATION 


A. OVERVIEW 

As the SRI experiment reported, the bomb moves at a high veloeity through the 
water and flow separation ereates a eavity of air around the body. That eavity, sometimes 
ealled cavitation, then remains in the water long after the bomb has passed and causes 
two areas of concern. First, will the trajectory remain stable, or will it tumble inside its 
own air cavity? Second, when the bomb does hit the cavitation wall, will the tail fins 
break? In addition, what happens to the trajectory after the fins break (Boeing 2009c)? 

Currently, three tests are available from which to draw data, as seen in Tables 2 
through 4. All three prove stability of the weapon to a certain depth in the water column, 
regardless of tail or fin separation. They also show that the tail fins most likely do (and 
the tail section possibly does) not remain intact during the full descent to the pond’s 
bottom, regardless of the impact AOA. 

According to the Boeing report (2009c), they found tail fin debris at the bottom of 
both UTT ponds. Though no underwater camera was available for UTT-1, once they 
drained the pond there was evidence that at least two tail fins came off the weapons. 
During UTT-2, camera images are only available for the first drop. They strongly 
suggest the first weapon’s tail impacted the cavity wall when the tail was about 1.5 m 
below the surface (the nose was around 5.5 m in depth), starting the process of breaking 
pieces off of the tail section, ultimately separating all four tail fins from the body. 

There is evidence that the second bomb also lost its fins. The SRI experiment 
shows that a bomb without tail fins is stable, but has significant curvature during its 
trajectory through the water. In the case of UTT-2(2), the bomb penetrated the bottom at 
such a shallow angle that it was able to burrow under the pond liner, climb the North face 
of the pond wall for a distance, re-enter the water traveling upward, and subsequently get 
airborne again (Boeing 2009c). This is a strong indicator that the bomb’s tail fins broke 
off far enough above the bottom of the pond to allow it to turn in the water. If this 
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shallow angle hypothesis is true, how is it that there is only a small horizontal deviation 
between the surfaee impaet point and the bottom impaet point? The following seetions 
will explain. 

B, MODEL VERIFICATION 

1. Lifting Surface Area Estimate 

STRIKE35 was ereated and tested on 1/12*'' seale replieas of the MK-84, so we 
must verify if the physics hold true for the full scale GBU-31, as well. As previously 
stated, there is only one model run that has the possibility of having maintained all four 
of its tail fins at a measurable point. LFFD-1 shows the best possibility of validating 
STRIKE35 because it traversed either without broken tail fins or the least distance with 
broken tail fins. An unknown, at the time of this validation, was the surface area of the 
lifting surfaces (tail fin and aero-surfaces or strake) (Figure 28). Known information 
about this drop is in Table 3. Another known fact is that STR1KE35 places the lifting 
surfaces solely at the tail section, so the validation needs to make up for this difference by 
changing the surface area of the tail fin in the model. 



Figure 28. GBU-31(V)2/B guidance set (U) [From Naval Air Weapons Center, Weapons 

Division 2002] 
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Because the data captured in LFFD-1 showed that the detonation occurred 12 
msec after water impact, so all lifting area comparisons were at 12 msec. Adjustments to 
the size of the fin surface area enabled an evaluation between the observed detonation 
horizontal deviation and the STRIKE35’s position. After changing the surface area of the 
fin and running the model, the nose and center of mass points helped solve for the 
equation of the body axis line. The known nose depth at detonation was plugged into the 
equation to solve for the horizontal displacement of the STRIKE35 weapon’s nose. The 
most accurate surface area, 0.68 m , brought the model within 6.8% of LFFD-l’s 
observed horizontal deviation, which is within the measurement uncertainty (Figure 29). 
This is much larger than the actual area of a tail fin pair (two tail fins in the same plane) 
and the aero-surfaces (strakes) combined (later found to be about 0.264 m ), but as a 
correction factor it works well as seen in the forthcoming validation. 

Fin Area Error 



0 2 03 04 05 06 0." Oi 

Fm .Aiea (in'^2) 

2 

Figure 29. The most accurate depiction of LFFD-1; single fin surface area of 0.68 m 

Calculating the accuracy for each “fin surface area” type is as follows. The 

displaced location, meaning the horizontal distance from initial water impact to the 
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detonation point, of the LFFD-1 weapon was 0.183 m from the point of water impact. In 
order to compare the STRIKE35 run (fin area of 0.68 m ), we found the equation of the 
line between the nose coordinate (0.4089, -4.048) m and the center of mass (0.04197, - 
2.487) m to be y = -4.2542x - 2.3085 (Figure 30). Plugging the detonation depth (y = 
3.14 m) in to this equation gives the calculated horizontal displacement (x) of the nose as 
0.195 m. This case is the least distance (0.012 m) between the observed and STRIKE35 
horizontal displacement. The error is the absolute value of the difference between the 
observed and the modeled horizontal deviation from water impact, all over the observed 
horizontal deviation (6.8%). 


tsrw. 0.012s. speed 321 



Figure 30. LFFD-1 with 0.68 m^ single fin surface area 
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2, Verifying Tail Breakage Physics of the Model 


Utilizing this adjusted fin area model, we ran all four of the remaining eomplete 
data sets (Tables 2 and 4) assuming the tail fins remained intaet (Figures 31 through 34, 
respeetively). All of the STRIKE35 runs were within one foot of the observed horizontal 
displacement, so it would seem that the weapons might have retained all four tail fins. 
However, there is strong evidence to the contrary. After draining the UTT-1 pond, pieces 
of tail fin were at the bottom. We also kn ow that during UTT-2(1) some of the tail fins 
definitively did not remain on the weapon. The underwater video cameras captured tail 
fin pieces leaving the cavitation wall when the nose was at a depth of about 5.5 m. 
Lastly, we saw in the SRI experiment that fewer tail fins caused more curvature in the 
weapon’s underwater trajectory; due to both UTT-2 weapons penetrating the deep pond’s 
bottom at a significant angle from vertical, this also lends to the hypothesis that both 
weapons lost tail fins. 


lime: 0.022s, speed: 315.2m/s 
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STRIKE35 0.429 m 


Observed 0.664 m 


Difference 0.235 m 
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Figure 31. UTT-1 (2) STRIKE35 model run with four fins 
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Figure 32. UTT-1(3) STRIKE35 model run with four fins 

tme: 0.03Ss, speed: 298 2in/s 



Figure 33. UTT-2(1) STRIKE35 model run with four fins 
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t»ne: 0.037s, speed: 290.6m/s 



Figure 34. UTT-2(2) STRIKE35 model run with four fins 


In order to use STRIKE35 to evaluate if and how many tail fins were lost in each 
run, we first needed to validate the model physics for the fins breaking mid-trajectory. 
We modified STRIKE35 (Appendix) to be able to input structural failure depth and the 
ability to choose either two fins or all four fins breaking off When the model runs UTT- 
2(2) without tail fins due to a structural failure at water impact, Figure 35 shows the 
weapon tumbles around 5 m. The horizontal displacement of the model (3.8 m at 12.2 m 
depth) is far worse than the observed. When UTT-1(2) and UTT-1(3) are modeled with 
the tail departing at water entry, they show very similar tumbling; the weapon tumbles 
beyond horizontal at the shallow pond bottom, and the horizontal displacement of their 
centroids is about 1.5 m at 7.6 m depth. By analyzing many more scenarios of tail 
breakage, the decision was made not to allow for tail failure input into STRIKE35 at this 
time. The amended program also does not add the option to select tail failure due to there 
being little evidence in the observed data of this actually occurring mid-water column. 
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tune: 0.076s. spe«t 79.7m/s 



The fact that we know the depth that UTT-2(1) lost a tail fin allowed us to 
extrapolate the possible scenarios. Analyzing two fins and four fins breaking off when 
the nose reached 5.5 m depth allowed a comparison against the different structural 
damages, as well as an evaluation of the model against the experiment. The observed 
data had a displacement from the water impact point of about 0.716 m. Figure 36 shows 
the model’s representation of two fins departing the weapon when the nose reaches 5.5 
m, and the results are within 0.061 m. The more accurate model depiction had four fins 
breaking off and a difference in horizontal displacement of 0.005 m from water impact 
point (Figure 37). Both demonstrations are within the measurement error of 0.3 m. 
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tw»e: 0 036$. $pee<t 298.4m/s 



time 0 036s, speed 297 8m/s 
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c. 


ASSUMPTIONS 


When creating STRIKES 5, they made several assumptions in order to make the 
model linear. When comparing the model to the observations, small errors could occur 
due to the following reasons; how the tail fins fall off the weapon, the shape of the nose 
cone, the model smoothes the body of the bomb, and the modeling location of two tail 
fins. When using the amended STRIKE35 version and selecting any tail fins to come off 
during the trajectory through the water, the fins fall off together and within a millisecond. 
The nose cone in JABS and the SRI experiments was the Navy Ogive type, whereas the 
SOABWFI Flight Testing experiments use the MXU-735 nose cone (Figure 38). 
Because STRIKE35 used the SRI data set in its creation, the model also uses the pointed 
Ogive nose cone, instead of the blunt MXU-735. The same reasoning applies to why the 
model does not include suspension lugs or the guidance band when creating the structural 
outline of the bomb (Figure 39). In order to maintain simplistic modeling features, they 
assumed that the weapon is an axially symmetric, rigid body; therefore, the two tail fin 
configuration is in the same plane instead of being perpendicular like the SRI experiment. 




Na\y Ogive MXU-735 

Figure 38. Ogive used for JABS, SRI and STRIKE35, whereas SOABWFI flight test 

used MXU-735 
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Lastly, utilizing the maximum density of pure water, 1000 kg/m^ at 4° C, is 
standard for seientifie models. The UTT-1 report stated the pond’s water temperature 
was about 35° C and the salinity near zero in the brackish pond (Boeing 2009b). This 
would put the density at about 996 kg/m^ (pure water at 35° C). Of course, the density 
difference will change the trajectory of the weapon through the water, as density is 
directly proportional to the hydrodynamic forces. Comparing the same run using 
STRIKE35 with its current density of 1000 kg/m^ to the observed density of 996 kg/m^ 
shows little difference though (Table 5). The average density for surface seawater (in the 
mixed layer) is about 1027 kg/m . Table 5 shows that the change in density from 1000 

3 3 

kg/m to 1027 kg/m makes a bigger difference in horizontal displacement than the 
change from 1000 kg/m to 996 kg/m (0.034 meters vice 0.004 meters), but both remain 
within the error of measurement. One interesting fact the table shows is that an alteration 
in density, though linear in the hydrodynamic forces, is not linear in the trajectory (or 
horizontal deviation at a certain depth). Instead, it changes the way the weapon “flies” in 
the water, and clearly in the ocean densities, it is not a linear change. 
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Water Type 

Density (kg/m^) 

Displacement (m) 

freshwater at 95° F 

996 

0.414 

freshwater at 4° C 

1000 

0.418 

saltwater 

1024 

0.462 

saltwater 

1025 

0.444 

saltwater 

1026 

0.390 

saltwater 

1027 

0.452 

saltwater 

1028 

0.450 

saltwater 

1029 

0.443 

saltwater 

1030 

0.432 


Table 5. Four fin GBU-31 nose displacement at 10 m depth. Water impact speed and 

AOA: 381 m/s and 1° 


D. MODEL APPLICATION 

1. Finding the Tail Breakage Depth 

In the cases used to validate STRIKE35 [LFFD-1 and UTT-2(1)], the model 
seems to work well, within measurement error, for either four fins or two fins breaking 
off the weapon. The three remaining cases we have complete impact data for [UTT-1(2), 
UTT-1(3), and UTT-2(2)] also have evidence that fins broke off while descending to the 
pond’s bottom, though it is unknown where and how many. For this challenge, we can 
use the adapted STRIKE35 for fins breaking in the water column and some logic to 
hypothesize how many fins and at what depth they broke. For each case, we ran 
STRIKE35, inputting the initial water impact characteristics, and broke off either two or 
all four tail fins at different depths (Table 6). Using the break-off option, two things are 
evident for the UTT experiments: 1) both types of breakage give results within 
measurement error, and 2) four fin failure shows the closest “observed to model” 
distances. 
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Experiment 

Breakage 

Depth 

(m) 

Model 

Displacement 

(m) 

Observed to 

2 Fin Failure 
Distance (m) 

Observed to 

4 Fin Failure 
Distance (m) 

UTT-1 (2) 

NA 

0.43 

0.24 

0.24 

UTT-1 (2) 

■■ 

0.44 

0.24 

0.23 

UTT-1 (2) 

3.3 

0.50 

0.23 

0.16 

UTT-1 (2) 

2.4 

0.62 

0.23 

0.05 

UTT-1(3) 

NA 

0.48 

0.09 

0.09 

UTT-1(3) 

■■ 

0.48 

0.09 

0.08 

UTT-1(3) 

3.3 

0.56 

0.09 

0.01 

UTT-1(3) 

2.9 

0.57 

0.09 

0.00 

UTT-2(2) 

NA 

1.23 

0.80 

0.80 

UTT-2(2) 

■■ 

1.33 

0.84 

0.90 

UTT-2(2) 

3.3 

2.30 

1.23 

1.87 


Table 6. Tail fin breakage at depth and assoeiated displaeements 


With all four fins breaking off at the same time, the least amount of error in 
bottom horizontal displaeement between the model and the observed occurred before the 
tail section hit the water for the UTT-1 experiment model runs. In UTT-1(2) and UTT- 
1(3), the nose depths were 2.4 m and 2.9 m, respectively. This is a strong indicator that 
the tail fins did not survive very long after impact with the water. Because the tail section 
does not enter the water until the nose is about 3.3 m, we have included the STRIKE35 
runs where all four tail fins break when the nose reaches 3.3 m depth of UTT-1 (2) and 
UTT-1(3) as Figures 40 and 41, respectively. The hypothesis that the fins departed not 
long after water impact directly contradicts the JABS trajectory analysis that showed an 
initial AOA less than two degrees will maintain structurally integrity beyond 28 msec; 
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not only did these two observations hit the water surface with less than 2° AOA, but 
STRIKE35 shows their trajectory to the bottom taking about 23 msec. 


tme: 0.023s, speed; 313.7tn/s 



Figure 40. UTT-1(2) with all four tail fins breaking off as the nose reaches 3.3 m 
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bme: 0 023$, spe«t 311 4m^ 



Figure 41. UTT-1(3) with all four tail fins breaking off as the nose reaehes 3.3 m 

UTT-2(2) was diffieult to figure out using Table 4 data and the amended 
STRIKE35. Working baekwards to find the tail breakage depth was more ehallenging 
due to the way the bomb travels in the water after a higher (greater than about 2°) water 
impaet AO A. A differenee between UTT-2 and UTT-1 is the pond depth, whieh one 
might think would give the weapon more time and even greater horizontal deviation at 
the pond bottom. Nevertheless, as seen in Table 6, breaking all four tail fins off at 
shallower depths inside the water eolumn created continuously growing horizontal 
deviation. Even keeping the tail fins and breaking only two fins off did not reduce the 
error. 

Another difference between UTT-2(2) and the previous UTT experiment is that 
the weapon was at a shallow enough angle when penetrating the pond bottom that is was 
able to tunnel. The tunneling behavior matches the SRI report claim that a bomb without 
tail fins has greater water trajectory curvature, as long as the weapon had enough depth to 
turn. Figure 42 shows that when all four tail fins break half way down the water column 
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(5.5 m depth in a 12.2 m pond), the bomb still does not have enough time to begin 
turning prior to hitting the bottom. From this, we infer that in order to get the bottom 
impact angle shallow enough to allow for tunneling, the bomb has to lose its tail fins 
closer to the surface. The SRI claim and the previous examination of UTT-1 leads to the 
guess that all four tail fins break off when the tail fins enter the water (when the nose 
reaches a depth of 3.3 m), but the error here is actually the highest (Table 6). 



Figure 42. UTT-2(2) with all four tail fins breaking off as the nose reaches 5.5 m 

After reading the previous paragraphs, it is curious how UTT-2(2) remained 
within 0.43 m of its water impact point. Interestingly, through a thorough examination of 
all initial inputs into STRIKE35 (with the density at 1000 kg/m^), without tail fins and 
above an initial impact AOA of about 2°, the trajectory through the water is extremely 
sensitive to impact AOA. Depending on the impact AOA, either the bomb’s trajectory 
through the water will continue smoothly in a curved path (as seen in the SRI data) or it 
will progress one way to a certain depth and then curve the opposing direction for the 
remainder of its path. The difference of 0.1° impact AOA can change the trajectory, and 
another 0.1° can change it back again. This turn-back trajectory can even happen at very 
low impact AOA (around 1 °), and is affected by impact speed, as well as the density of 
the water. 
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Detecting this sensitivity allowed for a supposition about the observation to model 
anomaly. With the observed impact AOA interpretation being ±0.5°, it is feasible that 
the impact AOA for the second weapon was slightly different from what the UTT-2 
report stated. If we increase the adapted STRIKES 5 (with all four tail fins removed at 3.3 
m) AOA input by 0.015° to 2.285°, the path of the bomb curves one direction and then 
turns to curve back towards the water impact point (Figure 43). This is the closest impact 
AOA to the UTT-2 report that results in the pitch-back trajectory, and the error is 0.70 m. 
Nonetheless, the closest model run to the observed initial characteristics and bottom 
impact point was the STRIKE35 initial AOA input increase of 0.025° to 2.295° (Figure 
44). This error is 0.30 m, which is within measurement error. 


tme: 0.038$, speed 261.6m/s 



Figure 43. UTT-2(2) when AOA increased to 2.285° 
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line 0.037s. speed: 276 4m7s 



2, Finding the Impact AOA Limitations 

The Technology Transition Agreement (TTA) between the Office of Naval 
Research and the Navy states that the “trajectory deviation from the flight path angle 
(velocity vector) should not exceed 7 ft (2.1 m) (approximate) for water depths of 10-40 
ft (3-12.2 m)” (Humes 2007). Based on these limitations, this study terminates analyzing 
different initial inputs when the weapon continuously exceeds 2.1 m from the water 
impact point. The forthcoming analysis also includes airspeeds within a reasonable 
window for a standard JDAM drop (381 to 396.4 m/s). Knowing that STRIKE35 
correlates well with LFFD-1 and UTT-2(1), both evaluated at a standard water density of 

-5 

1000 kg/m , we assume that the model also correlates when the density matches the 
environment of the observed weapon. 

The first impact AOA and speed analysis is of the China Lake 12.2 m pond, with 

•3 

a density of 996 kg/m . The second analysis is a hypothetical 12.2 m ocean environment 
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with a density of 1027 kg/m^. These eases are identieally analyzed as the worst seenario, 
with all four tail fins separating from the weapon upon water entry (when the nose 
reaehes 3.3 m) and the horizontal displaeement is measured from the water impaet point 
to the tip of the weapon’s nose, not the point of the explosive. The last analysis is of the 
detonation position, in the same ocean environment, for each of the delay settings within 
the FMU-139 fuze. These settings are 0 msec (impact), 10 msec, 25 msec, and 60 msec, 
but this report only includes the last three. The evaluation of each delay setting will be 
the same worst case described above. 

a. China Lake, Worst Case Analysis (Density of996 kg/m ) 

In Figure 45, the plot shows the horizontal displacement from the water 
impact point (m) as a function of AO A (°) in 35°C freshwater. Each line color represents 
a different water impact speed of the weapon, in increments of about 1.5 m/s. In every 
case, the weapon’s fins immediately fall off when the nose reaches a depth of 3.3 m and 
the weapon travels to a depth of exactly 12.2 m. This figure pictorially explains why 
UTT-2(2) was such a difficult case to analyze. The weapon can take a single trajectory, 
like the SRI experiment (Figures 10 through 12) or the pitch back trajectory, as seen in 
UTT-2(2) (Figures 43 and 44). 


57 



China Lake, Worst Case Analysis 



-•—381 m/s 
-*—382.5 m/s 

384.0 m/s 

- 385.6 m/s 

-iK- 387.1 m/s 
— • — 388.6 m/s 
—I— 390.1 m/s 
- 391.7 m/s 

393.2 m/s 
- 394.7 m/s 

396.2 m/s 


Figure 45. Weapon displacement vs. AOA in 35°C freshwater 


As seen by the coincidence of most lines below an impact AOA of about 
2°, the impact speed does not affect a weapon’s ultimate horizontal displacement from 
water impact point very much. Flence, impact speed does not change the weapon’s 
horizontal displacement as much as impact AOA does. Figure 45 also shows that when 
the weapon pitches back, it will not exceed the TTA until a higher impact AOA. Lastly, 
between about 1.2° and 2.1°, the displacements appear to increase linearly with the 
impact AOA (Figure 46). 
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China Lake, Worst Case Analysis 

—^381 m/s 
382.5 m/s 
384.0 m/s 

- 385.6 m/s 

^ 387.1 m/s 
-»— 388.6 m/s 

390.1 m/s 
— 391.7 m/s 

393.2 m/s 
394.7 m/s 

396.2 m/s 
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Figure 46. Most impact speeds have almost linear horizontal displacement between about 

1.2° and 2.1° 

All of the impact speeds remain within the TTA allowable deviation up to 
an impact AOA of 2.1°. At or below 384 m/s and at 387.1 m/s, the weapon continues the 
first type trajectory, exceeding the TTA between 2.1° and 2.2°. Between impact AOA of 
2.2° and 2.4°, these impact airspeeds transition to pitch back trajectory, allowing them to 
stay within the TTA longer. The best case for remaining within the limits of the TTA for 
the greatest range of AOA is an impact airspeed of 385.6 m/s. Of course, this is right in 
the middle (about ±1.5 m/s) of the impact speeds with lesser AOA tolerance, but this 
airspeed appears to have the pitch back trajectory type until the AOA of 2.6°. The 
remainder of the airspeeds transition to the pitch back trajectory at 2.1°, allowing these 
impact speeds not to exceed the TTA limit until 2.3° impact AOA. 

b. Ocean, Worst Case (Density of 1027 kg/m ) 

Figure 47 is the displacement (m) from the water impact point as a 

function of AOA (°) in an average mixed layer density in a hypothetical ocean. Each line 
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color represents a different water impaet speed of the weapon, in inerements of about 1.5 
m/s. In eaeh ease, the weapon’s fins immediately fall off when the nose reaehes a depth 
of 3.3 m and the weapon travels to a depth of exaetly 12.2 m. This figure shows that the 
weapon travels differently with an inerease in density. 


Ocean, Worst Case Analysis 


—^381 m/s 
382.5 tiVs 
* 384.0 m/s 

- 385.6 m/s 

-^387.1 m/s 
-•-388.6 m/s 
—I— 390.1 m/s 
^^391 .7 m/s 

393.2 m/s 
-394.7 m/s 

396.2 m/s 

1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 

Impact AOA (degrees) 

Figure 47. Weapon displaeement vs. AOA in a hypothetieal oeean mixed layer 

Just like the lower density ease, the mixed layer has the same two types of 
trajeetories. One readily observed differenee between the pond and the oeean is in the 
low AO As; all airspeeds, exeept 381 m/s, start out traversing the water with a piteh baek 
trajeetory. Between the impaet AOA of about 1.1° and 1.9°, the weapon transitions to the 
single eurve trajeetory style and remains within the TTA margins. Similar to the lower 
density ease, above 2° all impaet speeds experience flip-flop between both types of path. 
The best ease for remaining within the limits of the TTA for the greatest range of AOA 
are impaet airspeeds of about 387.1 to 390.1 m/s whieh allow an impaet AOA of up to 
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between 2.4° and 2.5°. The upper and lower limit of our sample, 381 m and 396.2 m, 
both only allow up to 2.1° AOA before departing the TTA boundary, while the remaining 
impaet airspeed allow an AOA of around 2.3°. 

c. FMU-139 Fuze, Worst Case (Density of 1027 kg/m^) 

The SOABWFI Flight Tests used the 10 msec delay on the FMU-139 fuze 
for the 7.6 m pond demonstrations (LFFD-1 and UTT-1) and the 25 msec delay for the 
12.2 m pond demonstrations (UTT-2). The last delay, not used yet in the SOABWFI 
program, is the 60 msec delay. This fuze also has a ± 20% tolerance that can detonate the 
bomb within the time limits. As in the previous examinations, the worst case means that 
all the tail fins break off when the nose is at 3.3 m (tail impact). During this 
investigation, we annotate the deviation from the water impact point as a function of time 
instead of depth. The model is run at each delay setting, at its lower limit (delay - 20%), 
and at its upper limit (delay + 20%). The horizontal displacement at the actual delay 
setting is the solid line, whereas the top of the error bar depicts the upper time limit’s 
displacement and the bottom of the error bar represents the lower tolerance’s deviation 
from water impact point. Because we are only interested in the horizontal deviation from 
the water impact point, the displacement at the actual delay setting is always the absolute 
value. Once again, the range of impact speed is from 381 to 396.4 m/s; this time in about 
3 m/s increments because the previous analysis shows the weapon is not as sensitive to its 
initial speed. 

In the ocean, the 10 ± 2 msec delay (Figure 48) remains within 1.8 m from 
the water impact point, no matter the impact airspeed or AOA. Similar to the findings of 
the SRI experiment, the greatest impact speed had the greatest deviation and the lowest 
impact speed had the least deviation from the water impact location for most of the 
impact AO As. One observation is that the greater the impact AOA, the greater the upper 
(+ 20%) and lower (- 20%) variation. In addition, at lower impact AOA (below about 
8°), the lower variation is much larger than the upper variation. 


61 



Detonation Displacement using the FMU-139 Delay Setting of 

10 ms + 20% 
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Figure 48. Weapon displaeement vs. AOA for 10 msec ± 20% delay 


The 25 ± 5 msec delay (Figure 49) remains within 2.1 m from the DPI, no 
matter the observed airspeed, up to an impact AOA of about 4°. Though the lowest 
impact speed (381 m/s) has a greater deviation from DPI below 3°, the highest speed 
(396.2 m/s) always has the greatest variation and has the highest deviation from 3° on. 
Opposite the shorter delay setting, this time the upper variation (+ 20%) seems to effect 
the deviation from DPI more so than the lower (- 20%). One observation is that the lower 
limit and the delay setting are always within the TTA, but it is the upper limit that 
reduces the available impact AOA. 


62 



























Detonation Displacement using the FMU-139 Delay Setting of 

2 5 ms + 2 0 % 
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Figure 49. Weapon displaeement vs. AOA for 25 msee ± 20% delay 


The 60 ± 15 msec delay (Figure 50) is not part of the TTA, but the 
availability of the setting makes it an interesting case study. All of the impact airspeeds 
and AOA up to 3° remain within about 8 m of the water impact point. One observation is 
that the higher impact speeds, 390.1 to 396.2 m/s, have the greatest variation in their 
tolerances. Allowing the user to choose less than 3° impact AOA and below 390.1 m/s 
impact speed, would reduce the horizontal deviation to 7 m. 
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Detonation Displacement using the FMU-139 Delay Setting of 

60 ms + 20% 
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Figure 50. Weapon displaeement vs. AOA for 60 msec ± 20% delay 
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V. DISCUSSION 


A. STRIKE35 UTILITY 

All of the live and inert weapons launehed into the pond at China Lake provide 
the proof of eoneept for the JABS as an organie mine neutralizer. The weapon maintains 
stability in the water column and remains within 2.1 m of the water impact point. 
STRIKE35 is useful to the proof of concept in that, with known start and end-points, it 
describes the path of the weapon through the water. The model can clarify if and where 
the weapon lost its tail section or tail fins. It is also valuable in validating other methods 
of data collection. 

In this analysis, STRIKE35 identified, using the impact points for the water and 
pond bottom, that it is probable the UTT tail fins did not remain intact long after 
submersion. Using it to investigate UTT-2(2) also shed light on the sensitivity to impact 
AOA, and the lack of sensitivity to impact speed of the weapon, on the overall horizontal 
deviation from water impact point. No matter the speed, about 2.1° is the cutoff impact 
AOA for predicting how the bomb will travel in the pond. Beyond this AOA, the weapon 
could continue to travel in a simple curve or could transition to a more complex, pitch 
back curve. The ability to predict which path the weapon will take becomes almost 
impossible, as a change as small as 0.015° impact AOA can transition the weapon from 
one path to the other. 

The model also shows that density affects the type of trajectory, especially around 
1 ° AOA, but does not affect the magnitude of the horizontal deviation from water impact. 
This is useful when applying the lessons learned about the weapon in the pond to the 
operational environment. Using a hypothetical ocean mixed layer density, STRIKE35 
then solved for the limitations on impact speed and AOA of the weapon for the FMU-139 
fuze delay settings. We found the 10 mecs delay is unlimited between 381 and 396.4 m/s 
impact speed, as well as beyond a 10° impact AOA. The 25 msec delay, for the same 
speeds, is limited to below a 4° impact AOA, and the 60 msec delay is currently not 
useful (under the TTA restrictions). 
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B, FUTURE USES AND RESEARCH 

Once the TTA is accepted, STRIKE35 ean assist ehoosing the parameters for a 
eorreet launeh aeeeptability region (LAR) when planning the AMCM mission. The LAR 
is a region where a launched weapon can successfully reach the target. Beeause the 
JDAM is a glide weapon, the kinematie performanee (energy of the weapon during flight) 
is the limiting faetor in defining the LAR. By seleeting the impaet of the water as vertieal 
(90°), the LAR shrinks to the point of only being able to drop a few weapons. Using 
STRIKE35, we observed the impaet AOA eould be less than vertieal and remain within 
the TTA limitations. 

STRIKE35 showed the limitations of the 60 msec delay fuze setting for this 
weapon and the TTA. SRI is eurrently looking at flatter nose cone shapes that ereate a 
larger eavity, allowing the weapon to travel straight for a longer distanee. If their 
experiments prove the eoneept, another amendment to STRIKE35 eould be useful in 
evaluating the new nose eone on the weapon. 

Some areas not looked at in this analysis ineludes the struetural integrity of the 
entire tail seetion, as well as wave effeets on the weapon’s horizontal deviation and 
struetural integrity. Aeeording to Boeing, the eonneetor bolt on the tail assembly is 
weaker than the fin attaehment bolts, but these three experiments were not affeeted by 
this fragile area. 
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VI. CONCLUSIONS 


According to ADM Roughead’s CNO Guidance for 2007-2008, “Adversaries are 
unlikely to attempt conventional force-on-force conflict, and to the extent that maritime 
forces could be openly challenged, their plans will almost certainly rely on asymmetric 
attack and surprise.” We have seen this in Iraq and Afghanistan with the improvised 
explosive device threat. In this same regard, the CNO has made the mine threat a top 
priority in his tenure. The reasons are elear; to be able to accomplish an amphibious 
landing, keep sea lines of eommunieation and commerce open, deny the threat to a port 
disrupted by a mine laying campaign, and maintain freedom of maneuver and freedom of 
access. 

Organic mine warfare is the newest form of mine warfare. In order to bridge the 
gap between the aging AMCM and the latest technology, a TTA took place to allow for 
the use of the GBU-31 as a mine neutralizer. Resources are in high demand and short 
supply to be able to test all possible scenarios, so the STRIKE35 model can assist with 
the final stages of the TTA. 

This thesis validates STRIKE35, which is based on the 1/12*’’ seale repliea of the 
MK-84, for the GBU-31. After validation, an amendment to the model allowed us to 
work backwards to find the depth where the tail fins broke off for three other weapons. 
One of the trajectories showed the sensitivity that the JDAM has to impact AO A. The 
AO A limitations found for a 12.2 m freshwater pond, as well as a hypothetical ocean 
(with the density 1027 kg/m^) of the same depth, were below 2.1°. We found the impact 
speed and water density to have little bearing on the overall horizontal trajectory of the 
weapon from the water impaet point. Finally, we observed speeific fuze settings and 
their tolerances in the ocean environment, finding the 10 msec fuze to have no limitations 
and the 25 msec to be limited below 4°. 
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APPENDIX 


runMKSDtailbreak.m 

% this program has changed STRIKE35 to show what happens when the tail fins 
% break off inside the water column, after water entry. 

% Created by Chenwu Fan 

global g rw m mw mwfl J L Le d nu area yi areai voli veyi emr CGy yyy ReO nfin 
% set up constants 

dist=0.1; % the distanee for eaeh step 

ipt=round(2/dist);% how often to plot figure of bomb 
r2d=l80/pi; g=9.806; % radians to degrees; gravity; 
f2m=.3048; % foot to meter or ft/see to m/s 
nu=1.13e-6; ftsz=8; % dynamatie viseoeity; font size 
rw=1000; % density of water 

mfl=6.2425*0.454; volfl=mfl/2200;% mass of fin (kg); volume of the fin 
figure('unitsVinchesVposition',[2,2,8,6]);% opening figure 
% eharaeter array (no tail, no fins, 2 fins, 4 fins) 

% Assume Original Configuration is 4; 

tt=input('Enter Water Depth (ft). Fin Break-Off Depth (ft), MK84 Break-Off Configuration fin s : (0,2), Speed (ft/see). Impact Angle of 
Attaek (degrees) 

while(~isempty(tt))% as long as eharaeter array input, this loop will run, else it exits program 
nfin=4; b=0; % tail not broken 

elf;% elears eurrent figure 

tt=str2num(tt);% eonverts a eharaeter array of numbers to a numerie matrix 
dep=tt(l); fmbreak=tt(2); bnfm=tt(3); vs0=tt(4); EL0d=tt(5); % assigns names to inputs 

dep=-l*dep*f2m; linbreak=-l*finbreak*f2m; vs0=vs0*f2m; EE0=EE0d/r2d; % unit conversion from english to metric 
n=l; 

load MK84tailsize.dat;% outline of bomb shape 

yi=MK84tailsize(:,l); ri=MK84tailsize(:,2);% length from nose to tail; radius at respeetive points 
E=150.51*0.0254; % length converted to em 
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m=2076.64*0.454; CGy=63.12*0.0254;% mass converted to kg; eenter of gravity eonverted to em 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

[voli,areai,voyi,Jy,Jz,omd,ovd,yi,ri,yoi,omr]=getbombphysios(yi,ri,CGy,m);% program for how the bomb travels through the water 
vol=voli(end); mw=rw*vol; d=2*sqrt(vol/L/pi);% renames total volume; mw=water mass=density*volume; diameter of eylinder of 
equal volume 

J = Jz; area=areai(end);% renaming inertia and total area 

m=m+nlin*mfl; mwfl=nfin*volfl*rw; mw=mw+mwfl; % if fins then add the mass and volume to bomb, and mass of water 
displaeed 

LH=yoi(l); % first y variable based on mass center and head forward 

LT=L-LH; Lo=vcyi(end)-omd;% length minus previous; volume center locations along the bomb-center of mass 
% finding points for tail dimensions 

yl=yoi(end)+0.75*0.0254; y2=yl; y3=yoi(end)+(0.75+9.44)*0.0254;% distanee from head eonverted to em 

rl=interpl(yoi(end-l;end),ri(end-l:end),yl);% finds the radius at the end point 

r2=25.31/2*0.0254; r3=r2;% eonverts the radius to em 

yy=yoi(end-l;end); rr=ri(end-l;end);%defmes the eurrent last two points 

dy=diff(yy); dr=diff(rr);%find the slope between eurrent last two points 

y4=(dr*yy(l)-l-dy*(y3+r3-rr(l)))/(dy+dr);% defines a new last point 

r4=y3+r3-y4; % defines a new last radius 

yyy=(yl+y2-i-y3+y4)/4;% average 

yfil=[yl,y2,y3,y4,yoi(end)]; rfil=[rl,r2,r3,r4,ri(end)];%fills in the length and radii 
yyi=[yoi,yfil,fliplr([yoi,yfil])]; rri=[ri,rfil,-fliplr([ri,rfil])];%ereates bottom outline of the bomb 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

axes('position',[0.1,0.1,0.85,0.85],'ylim',[dep-0.5,4.5],'boxVonVfontsize',ftsz);% set plot oharaeteristies 

hold on; axisCequaf); 

plot([-l,10],[0,0],'b-'); 

plot([0,0],[dep,0], 

ylabel('Y(m)'); xlabel('X(m)'); 

elear XYZC EL AG VS;% clears center eoordinate, lift direetion, angle of velocity, bomb speed 
OMG0=zeros(l,3); % zeros plaeeholder 

psiv=pi/2; psib=pi/2;%bomb veloeity and bomb position azimuth angles set to 90 degrees 
EL(l)=EL0-pi/2;% rotates bomb position elevation angle 90 degrees 
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VS(l)=vsO; AG(l)=-pi/2;% initial bomb speed and bomb veloeity elevation angle 
XYZC(:,l)=[0;0;CGy];% replaeing first column of center coordinate matrix 
[myi,mri]=rotobj'(yyi,rri,EL(n));% the object (yyi,rri) rotation arf (radiau) around (0,0) 
hh2=patch(XYZC(2,n)+myi,XYZC(3,n)+mri,'y'); % colors the bomb yellow 
plot(XYZC(2,n),XYZC(3,n),'r.');% plots the center of mass in red 
Lw=0; time(n)=0; ELAG=0;% initial 

omg3D=[le-10;0;0]; kwet=0; % initial bomb rotation velocity; kwet 

while(XYZC(3,n)+min(mri)>dep & (kwet==0 | XYZC(3,n)+max(mri)<0))%while the bomb has not reached depth 
dt=min( 1 ,dist/(VS(n)+eps));%time step 
nl=n; n=n+l; time=time+dt;%initial time 

[VS(n),AG(n),EE(n),XYZC(:,n),omg3D,psiv,psib]=bomblstep3D(dt,VS(nl),AG(nl),EE(nl),XYZC(:,nl),omg3D,psiv,psib); 

[myi,mri]=rotobj(yyi,rri,EE(n));% the object (yyi,rri) rotation arf (radiau) around (0,0) 

if(psib<0), myi=-myi; end% if bomb position azimuth angle negative then places bomb on other side 

if(XYZC(3,n)<-2.5), kwet=l; end% if depth less then -2.5 then wet 

if(ipt*floor((nl-l)/ipt)~=nl-l), set(hh2,'visible','off); end 

hh2=patch(XYZC(2,n)+myi,XYZC(3,n)+mri,'y');% colors the bomb yellow 

plot(XYZC(2,n),XYZC(3,n),'r.');% plots the center of mass in red 

title(['time: ',sprintf('%7.3f,time(end)),'s, speed; ',sprintf('%6.1f,VS(n)),'m/s']); 

drawnow; 

if(b==0 & XYZC(3,n)<=finbreak) 
b=l; nfin=bnfin; 

text(XYZC(2,n),XYZC(3,n),'-','color','k','horizontalalignmenf,'center'); 

if(nfin==0)% no fins 

m=2051.67*0.454; CGy=62.11*0.0254;% mass converted to kg; center of gravity converted to cm 
elseif(nfin==2)% two fins 

m=2064.15 *0.454; CGy=62.60*0.0254;% mass converted to kg; center of gravity converted to cm 
% elseif(nfin==4)% four fins 

% m=2076.64*0.454; CGy=63.12*0.0254;% mass converted to kg; center of gravity converted to cm 

end 

m=m+nfin*mfl; mwfl=nfin*volfl*rw; mw=mw+mwfl; % if fins then add the mass and volume to bomb, and mass of water 
displaced 
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if(nlin==0) 

yyi=[yoi,fliplr(yoi)]; rri=[ri,-fliplr(ri)];% update bomb shape 
end 
end; 
end 

tt=input('Enter Water Depth (ft), Fin Break-Off Depth (ft), MK84 Break-Off Configuration fins; (0,2), Speed (ft/sec), fmpaet Angle 
of Attaek (degrees) 
end 

funetion [ee,ev,ealpha,el,ed,evgama,evpsi]=setunitveetor(gama,beta,psiv,psib); 

% function [ee,ev,ealpha,el,ed,evgama,evpsi]=setunitvector(gama,beta,psiv,psib); 

% is used calculate unit vector 
% input; 

% gama; bomb velocity elevation angle. 

% beta; bomb position elevation angle. 

% psiv; bomb velocity azimuth angle 
% psib; bomb position azimuth angle 
% output; 

% ee; bomb axis direction vector 

% ev; bomb velocity direction vector 

% ealpha; attach angle (Hydro torque) direction vector 

% el; lift force direction vector 

% ed; drag force direction vector 

% evgama, evpsi; d(ev)/dt two compound vector. 

% 

% ChenwuFan, NOAP, NPS 

ee=[cos(beta).*cos(psib); cos(beta).*sin(psib); sin(beta)]; % bomb direction 
ev=[cos(gama).*cos(psiv); cos(gama).*sin(psiv); sin(gama)]; % velocity direction 
ealpha=cross(ev,ee); 

ealpha=ealpha./(sqrt(sum(ealpha.^2))+eps); % torque direction 

el=cross(ealpha,ev); % lift direction 

ed=-ev; % drag direction 
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evgama=[-sin(gama).*cos(psiv);-sin(gama).*sin(psiv);cos(gama)]; % d(ev)/dt gama 
evpsi=[-sin(psiv); cos(psiv); zeros(size(gama))]; % d(ev)/dt psiv 


function [Cd,Cl]=syniNACAOO 15coef(ang,dor); 

% function [Cd,Cl]=symNACA0015coef(ang,dor); 

% dor; 1-units degree, 0-units radiu (default). 

% ChenwuFan, NOAP, NPS 
if(~exist('dor') | -dor), ang=ang* 180/pi; end 

agel=[0,0;l,0.045;2.5,0.23;4.7,0.5;6,0.7;9,0.9;10.05,0.95;12.5,1.055;14,0.95;15,0.62; 

16,0.6;17,0.59;18,0.585;20,0.6;22,0.61;25,0.68;30,0.88;35,0.99;40,1.02;45,1.03; 

50,1;55,0.94;60,0.86;65,0.75;70,0.6;75,0.47;80,0.35;85,0.21;90,0.09;93,0; 


100,-0.18;105,-0.3;110,-0.42;115,-0.58;120,-0.7;125,-0.785;130,-0.88;135,-0.98; 

140,-0.98;145,-0.9;150,-0.78;152.5,-0.7;155,-0.66;160,-0.62;165,-0.695;170,-0.89; 

175,-0.67;180,0]; 

aged=[0,0;10,0.01;15,0.08;20,0.25;25,0.4;30,0.59;35,0.71;40,0.9;45,l.l;50,1.24; 
55,1.4;60,1.5;65,1.6;70,1.67;75,1.72; 80,1.76;85,1.8;90,1.81;95,1.78;100,1.73; 
105,1.68; 110,1.62;! 15,1.55;120,1.5;125,1.4;130,1.27;135,1.1;140,0.95;145,0.76; 
150,0.6;155,0.4;160,0.31;165,0.22;170,0.12;175,0.03;180,0.01]; 

Cd=interp 1 (aged(:, 1 ),aged(; ,2),ang); 

Cl=interp 1 (agel(:, 1 ),agel(: ,2),ang); 


funetion F3d = onefoilforee(vf3d,ee); 

% funetion F3d = onefoilforoe(vf3d,ee); 
% vf3d: fin eenter veloeity (3D) 

% ee; bomb direetion (unit veetor) 

% F3d; 3-D foree on one fin 
% C. Fan, NOAP NPS 9/21/2009 
rw=1000; finarea=0.68^2; 
vs=sqrt(sum(vf3d.^2)); 
if(vs<eps) 

F3d=zeros(3,l); 

return; 
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end 

ev=vf3d/vs; 

alpha=acos(sum(ev.*ee)); 

ealpha=cross(ev,ee); 

ealpha=ealpha./(sqrt(sum(ealpha.^2))+eps); % torque direetion 
el=cross(ealpha,ev); % lift direetion 

ed=-ev; % drag direetion 

[Cd,Cl]=symNACA0015ooef(alpha); 

Fd=0.5 *Cd*rw*linarea. *vs A2; 

Fl=0.5*Cl*rw*linarea.*vsA2; 

F3d=Fd*ed+Fl*el*0.6366; _ 

function [OMG,IntOMG] = OMGsolve(lm,B,OMGO,dt); 

% function [OMG,IntOMG] = OMGsolve(lm,BB,OMGO,dt); 

% Solve diagonal ODE as dOMG/dt = lm*OMG+ B, with the initial 
% condition OMGO at time step dt. 

% 

% ChenwuFan, NOAP, NPS 4/10/2007 
if(lm<-le-8) 

tenninal=-B/lm; c=OMGO-terminal; 

OMG=c*exp(lm*dt)+terminal; 

IntOMG=terminal*dt+c/lm*(exp(lm*dt)-l); 

else 

OMG=OMGO+(B+lm*OMGO)*dt; 

IntOMG=OMGO*dt+(B+lm*OMGO)*dG2/2; 

end 

end 

function Cm = MK84cm(elag,Re,del); 

% ChenwuFan, NOAP, NPS 

refRe=l .8e7; % Critical Reynolds Number 

re=Re/refRe; 

Cm=NaN*zeros(size(re)); 
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il=find(elag<=pi/2); 

Cm(il)=0.07*sin(2*elag(il))./re(il).^0.2; 

i2=find(elag>pi/2); 

Cm(i2)=0.02 *sin(2 *elag(i2)). *re(i2) .^0.5; 
function cl=MK84cl(elag,Re,del); 

% ChenwuFan, NOAP, NPS 
% elag=abs(elag); 

refRe=l .8e7; % Critical Reynolds Number 

re=Re/refRe; 

cl=NaN*ones(size(Re)); 

il=find(elag<=pi/2); 

afl=pi*(2*elag(il)/pi).^l .8; 

cl(il)=0.35*sin(afl).*re(il).^0.2; 

i2=find(elag>pi/2); 

af2=2*pi*(2*elag(i2)/pi-l).^0.7; 

cl(i2)=0.1*sin(af2)-0.015*del(i2).*re(i2).^2.*sin(af2.^0.85); 
function cd = MK84cd(elag,Re,del); 

% ChenwuFan, NOAP, NPS 

% according to the scale model experiment. Re > 3x10^5 for bomb velocities 

% > 7.80288 m/sec and drag eoeffieient < 0.2 

refRe=l .8e7; % Critical Reynolds Number 

re=Re/refRe; 

af=pi-2*elag; 

% aff=sign(af).*sqrt(pC2-(pi-abs(af)).^2); 

% ed=0.02+0.25*exp(-(elag-pi/2).^2*2).*re.^0.2+0.008*sin(aff).*del; 

N=2.2; 

aff=sign(af).*(pCN-(pi-abs(af)).^N).^(l/N); 

cd=0.02+0.35*exp(-(elag-pi/2).^2*2).*re.^0.2+0.008*sin(aff).*del; 

function [voli,areai,vcxi,Jx,Jy,cmd,cvd,xi,ri,xoi,cmr]=getbombphysics(xi,ri,CGx,m,denrat); 

% function [voli,areai,vcxi,Jx,Jy,cmd,cvd,xi,ri,xoi,cfr,cmr]=getbombphysics(xi,ri,CGx,m,denrat); 
% is used to ealeulate bomb physics parameters. 
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% input: 

% xi,ri: the Bomb Radius data set (xi: distanee from head) 

% CGx: eenter of Gravity (or mass) from front point (m) 

% m: mass of bomb (kg) 

% denrat: density ratio distribution as:[x,x2,x3,...;fao,fao2,fao3,...], default: [0,L;1,1] 

% output: 

% voli: volume data set 

% area: eonjeetion area data set 

% vexi: volume eenter data set (distanee from head) 

% Jx,Jy: Inertia around x, y axis 

% emd: eenter of mass distanee, from volume eenter to mass eenter 
% (>0 Me ahead of volume eenter) 

% evd: from volume to middle point. (>0 middle ahead of Ve) 

% xi,ri: similar as input but add mass eenter point 
% xoi: X variable based on mass eenter and head forward 

% emr: the anti-rotation torque eoeffieient Mr = -0.5*mw*C*|omg|*omg*L^2*emr 
% ChenwuFan, NOAP, NPS 

if(~exist('denraf)), denrat=[0,xi(end);l,l]; end % if the density ratio not input, then given a value of 1 
if(size(xi,l)>l), xi=xi'; ri=ri'; end % if bomb has a length then ehange row veetor into eolumn veetor 

iil=lind(xi<=CGx); ii2=lind(xi>=CGx); % finds and names the length that's less than or equal to and greater than or equal to 

the eenter of gravity 

if(iil(end)~=ii2(l)) % if they are not equal, then the radius must be found at the eenter of gravity 

ro=interpl(xi,ri,CGx); % the radius is interpolated at the eenter of gravity 

ri=[ri(iil),ro,ri(ii2)]; % the interpolated radius is then added to the radius matrix 

xi=[xi(iil),CGx,xi(ii2)]; % as well, the eenter of gravity is added to the length matrix 

end 

dxi=diff(xi); % finds the distanee between eaeh point of length matrix 

areai=(ri(l :end-l)+ri(2:end)).*dxi; % begins finding the 2D area by integrating the radius aeross the length of the bomb 

areai=[0,areai]; areai=oumsum(areai); % adds a zero/start point to the matrix; gives the total area in sliees from beginning to 

end of bomb 

rl l=ri(l:end-l).^2; rl2=ri(l:end-l).*ri(2:end); r22=ri(2:end).^2; % radii squared 
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dri=diff(ri); % finds the distance between each radius 

voli=pi*(rl l+rl2+r22)/3.*dxi; % frustum volumes (cylindar volume = pi*r^2) 

vci=(rl l+2*rl2+3*r22)./(4*(rl l+rl2+r22)).*dxi; % volume center of each slice 

vcxi=vci+xi(l;end-l); % volume center locations along the bomb 

faci=interpl(denrat(l,:),denrat(2,:),vcxi); % volume center locations with density ratio factored in 

Vco=sum(voli.*vcxi)/sum(voli); % volume center (from nose of bomb) 

cmd=Vco-CGx; cvd=Vco-xi(end)/2; % center of mass=volume center-center of gravity;from volume to middle point 

xoi=CGx-xi; % coordinate system moved to center of gravity (x variable based on mass center and head 

forward) 

[mmm,io]=min(abs(xoi)); % finds the [value,location] of the minimum in the coordinate system 

idxf=[l ;io-l]; idxb=[io;length(vci)]; % creates front and back shapes 

volf=sum(voli(idxf)); volb=sum(voli(idxb)); % creates front and back volumes 

volxf=sum(voli(idxf).*faci(idxf).*(vcxi(idxf)-CGx));% front volume center locations with density ratio factored in 

volxb=sum(vob(idxb).*faci(idxb).*(vcxi(idxb)-CGx));% back volume center locations with density ratio factored in 

rhof=m*volxb/(volf'volxb-volb*volxf); % density of front 

rhob=m*volxf/(volb*volxf-volxb*volf); % density of back 

rhoi=rhof‘'faci; % density of back with density ratio factored in 

rhoi(idxb)=rhob*faci(idxb); 

vcxi=cumsum(vcxi.*voli)./cumsum(voli); vcxi=[0,vcxi];% the average cumulative sum of volume center locations along the bomb 

voli=[0,voli]; voli=cumsum(voli); % the cumulative sum of frustum volumes along the bomb 

Jxi=(rl 1.^2+rl2.*rl 1+rl I.*r22+rl2.*r22+r22.^2).*dxi.*rhoi*pi/10; % polar moment of inertia in x of sub-frustum 

Jx=sum(Jxi); % x-total polar moment of inertia 

drdx=dri./diff(xoi); xl=xoi(l;end-l); x2=xoi(2;end); ril=ri(l;end-l); 

ridr=ril-drdx.*xl; 

JJi=ridr.M.*(xl-x2)/4+ridr.^3.*drdx.*(xl.^2-x2.^2)/2+ridr.^2.*(l/3+drdx.^2/2).*(xl.^3-x2.^3)+... 

ridr.*drdx.*(l/2+drdx.^2/4).*(xl.M-x2.M)+drdx.^2.*(l/5+drdx.^2/20).*(xl.^5-x2.^5);% polar moment of inertia in y of sub¬ 
frustum 

Jy=pi*sum(rhoi.*JJi); % y-total polar moment of inertia 

% The anti rotation torque caused by rotation as Mr = -0.5*mw*C*|omg|*omg*L^2*cmr 

xm=(xl+x2)/2; 
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iil=find(xm>0); 

cmrl=ridr(iil).*(xl(iil).M-x2(iil).M)/4+drdx(iil).*(xl(iil).^5-x2(iil).^5)/5; 

ii2=lind(xm<0); 

cmr2=ridr(ii2).*(xl(ii2).M-x2(ii2).M)/4+drdx(ii2).*(xl(ii2).^5-x2(ii2).^5)/5; 
cmr=2*[cmrl,-cmr2]; cmr=cumsum(cmr); 
cmr=[0,cmr]/(voli(end)*xi(end)^2); 


function cd2 = CdR28(Re); 

% function cd2 = CdR28(Re); 

% Impact 28 Cross direction drag coefficient. 

% 

% C. Fan NOAP NPS 9/10/2003 
if(Re>=350000) 
cd2= 1/(641550/Re + 1.5); 
elseif Re>= 150000 
ed2= 1.875 -.0000045 * Re; 
elseif Re>= 12000 
cd2 = 1.2; % aspect ratio? 

elseif Re>= 2000 
ed2 = .84 + .00003 * Re; 
elseif Re >= 180 

cd2 = .86 + 89 / Re; % .855 +89/Re 
elseif Re >=12 
cd2= 1.26+ 16/Re; 
elseif(Re>le-10) 

ed2 = 1.8 + 8 / Re; % 1.9276 + 8/Re 
else 
cd2=0; 
end 

function [VS,GAMA,BETA,XYZ,OMG3D,psiv,psib]=bomblstep3D(dt,vs,gama,beta,xyz,omg3D,psiv,psib); 
% beta(el); bomb position elevation angle. 
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% gama(ag): bomb velocity elevation angle. 

% psiv: bomb velocity azimuth angle 
% psib; bomb position azimuth angle 
% omg3D(del); bomb rotation velocity 
% ChenwuFan, NOAP, NPS 

global g rw m mw mwfl J L Lc d nu area yi areai voli vcyi cmr CGy yyy ... 

ReO nfin 
fff=1.85; 

if(nfm==0), fff=0.7; end 
% get 3D unit vectors 

[ee,ev,ealpha,el,ed,evgama,evpsi]=setunitvector(gama,beta,psiv,psib); 

% calculate attached angle 
alpha=acos(sum(ev.*ee)); 

Re=vs*d/nu; % initial reynolds number = initial velocity/(diameter of cylinder of equal volume*dynamic viscocity) 

Lw=-xyz(3)+L/2; Lw=mm(Lw,L); Lw=max(0,Lw); % Under water bomb length 

omg=sum(omg3D.*ealpha); % angle velocity attached angle component 

% calculate hydro dynamic coefficients; only drag coefficient effected by 

% impact factor 

Acd=MK84cd(alpha,Re,omg); 

Acl=MK84cl(alpha,Re,omg); 

Acm=MK84cm(alpha,Re,omg); 
factimp=mm(3,L/(Lw+eps)); % impact factor 
Acd=Acd*factimp; 

% As hydro force only effect the underwater part, must calculate the 

% underwater project area(areaa), volume(voll), volume center(vcyy), and 

% mass center (cmrr). 

areaa=interp 1 (yi,areai,Lw); 

volUinterp 1 (yi,voli,Lw); 

vcyy=mterp 1 (yi,vcyi,Lw); 

cmrr=interp 1 (yi,cmr,Lw); 

Fbuo=voll*rw*[0;0;g]; % 3D buoyancy force 
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rav2=rw*area*vs^2; 

Fd=0.5*Acd*rav2; Fl=0.5*Acl*rav2; 

F3dd=Fd*ed+Fl*el; % 3D hydro dynamic force 

F3d=F3dd+m*[0;0;-g]+Fbuo; % total 3D force (hydro, gravity andbuo forees) 
Mt3d=0.5*Acm*rav2*L*ealpha; % Hydro dynamic torque 
MOd=cross((CGy-vcyy)*ee,F3dd); % 3D hydro dynamic force torque 
M3d=Mt3d+Mf3d; 
if(nfm>0) 

if(xyz(3)<Lc-L) 

Mfbuo=cross(yyy*ee,mwfl*[0;0;g]); % fm buoyancy torque 
Fbuo=Fbuo+mwfl*[0;0;g]; % add fm buoyaney 
F3d=F3d+mwfl*[0;0;g]; % add fm buoyancy 
Ff3D=onefoilforce(vs*ev+cross(omg3D,yyy*ee),ee)*nfm; 

F3d=F3d+FBD; % add fm 3D forces 
Mf3df=cross(yyy*ee,FOD); 

M3d=M3d+(Mf3df+Mfbuo); % add fm 3D torque 
end 
end 

M3d=M3d+eross((CGy-veyy)*ee,Fbuo); % add torque eaused by Buoy 

% As 3-D, M3d need seperate into OMG and (OMG x ee) direction componemts. 

OMG=sqrt(sum(omg3DA2)); 

if(OMG<le-20), eomg=[l;0;0]; else, eomg=omg3D/OMG; end % direetion of OMG 
eomgME=oross(eomg,ee); 

eomgME=eomgME/sqrt(sum(eomgMEA2)); % direction of OMG x ee 
Vax=sum(ev.*ee)*vs; 

Vr3D=vs*ev-Vax*ee; 

Vr=sqrt(sum(Vr3DA2)); 

MrJ=-CdR28(Vr*d/nu)*rw*d*E^3*Vr/12/J/fff; % anti rotation 
% Caleulate angle inerease in OMG direetion 
M J=sum(M3 d. * eomg)/ J; 

[OMG,IntOMG] = OMGsolve(MrJ,MJ,OMG,dt); 
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IntOMG3D=IntOMG*eomg; OMG3D=OMG*eomg; 

% calculate angle increase in normal OMG direction (OMG x ee) 

MMEJ=sum(M3 d. * eomgME)/J; incangME3 D=0.5 *MMEJ * dG2 * eomgME; 
totlincang3 D=IntOMG3 D+incangME3 D; 
ebeta=[sin(psib);cos(psib);0]; 
dbeta=sum(totlincang3D. *ebeta); 

BETA=beta+dbeta; psib=psib+totlincang3D(3); 

As=sum(E3d.*ev)/m; % speed aeeelerate 
Dgama=sum(E3d.*evgama)/(m*vs); % d(gama)/dt 
Dpsiv=sum(E3d.*evpsi)/(m*vs); % d(psi)/dt 
agg=gama+Dgama*dt/2; % middle step gama 
psivm=psiv+Dpsiv*dt/2; % middle step psiv 

evm=[cos(agg).*cos(psivm); cos(agg).*sm(psivm); sm(agg)];% middle step ev 
% update XYZ, gama & psiv 
X YZ=xyz+( VS+ As * dt/2) * dt* e vm; 

VS=vs+As*dt; GAMA=gama+Dgama*dt; 
psiv=psiv+Dpsiv*dt; 

% limit BETA & GAMA value with in (-pi,pi] 
while(BETA>pi), BETA=BETA-2*pi; end 
while(BETA<=-pi), BETA=BETA-l-2*pi; end 
while(GAMA>pi), GAMA=GAMA-2*pi; end 
while(GAMA<=-pi), GAMA=GAMA-l-2*pi; end 
% update BETA and GAMA value within [-pi/2,pi/2] 
if(BETA>pi/2), BETA=pi-BETA; psib=psib+pi; end 
if(BETA<-pi/2), BETA=-pi-BETA; psib=psib+pi; end 
while(psib>pi), psib=psib-2*pi; end 
while(psib<=-pi), psib=psib+2*pi; end 
if(GAMA>pi/2), GAMA=pi-GAMA; psiv=psiv+pi; end 
if(GAMA<-pi/2), GAMA=-pi-GAMA; psiv=psiv+pi; end 
while(psiv>pi), psiv=psiv-2*pi; end 
while(psiv<=-pi), psiv=psiv+2*pi; end 
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